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FOREWORD 


The PLANS system of programs for the nonlinear analysis of 
structures was developed within -the Grumman Research Department 
with the cooperation, guidance, and partial support of the NASA 
Langley Research Center. 

The programs are an outgrowth of the work reported in a 
series of NASA Contractor Reports: CR-803, CR'1649, CR-2351, and 

CR-2568. The last-named document is the theoretical companion 
volume to this manual. 

The principal developers of PLANS are Drs. A. Pifko, H. Armen Jr 
H. Levine, and A. Levy. The successful development of a system as 
broad in scope and complexity as PLANS requires the efforts of many 
individuals. The principal developers gratefully acknowledge the 
contributions of J. S. Miller for his efforts associated with the 
initial programming. Special thanks also go- to Ms. P. Ogilvie, 

P. Zirk, and E, Yander for their diligence and dedication to the 
programming effort. 
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SECTION 1 


INTRODUCTION AND INPUT OVERVIEW 



INTRODUCTION AND INPUT OVERVIEW 


The PLANS system, rather than being one comprehensive computer 
program, is a collection of finite element programs used for the 
nonlinear analysis of structures. This collection of programs 
evolved and is based on the organizational philosophy in which 
classes of analyses are treated individually based on the physical 
problem class to be analyzed. On the basis of this concept, each 
of the independent finite element computer programs of PLANS, with 
an associated element library, can be individually loaded and used 
to solve the problem class of interest. A number of programs have 
been developed for material nonlinear behavior alone and f.or com- 
bined geometric and material nonlinear behavior. Table 1 summarizes 
the usage, capabilities, and element libraries of the current pro- 
grams of the PLANS system. These include: 

* BEND for the plastic analysis of 

built-up structures where bending 
and membrane effects are sig- 
nificant 

* HEX for problems requiring a three 

dimensional elastic-plastic 
analysis 

a REVBY for the plastic analysis of 

bodies of revolution 

6 OUT-OF-PLANE MG for the material and geometric 

nonlinear analysis of built-up 
structures 

Supplementing these is a SATELLITE program for data debugging 
and plotting of input geometries. 


1-1 



1-2 


TABIE 1 OPERATIONAL MODULE S OF PLANS - CURRENT STATUS (5/77) 
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no 
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Ho 

Yes 
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Nonlinearity (Sc3U 
strain-large deflect.) 

• Bear (9 cross- 
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« Stringer 

• Triangle Merbrur.e 
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9 Hiersal (E*. 

Scaa„) 
o Cyclic 

• Non - Proper tionoi 

a As In 
BSNC 

Yec 

No 

Yea 

ft* 

(see 
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In addition, a number of spin-off special purpose modules 
have been developed: FAST for fracture analysis, AXSHEL for the 

nonlinear analysis of shells of revolution, and BENST for the 
buckling analysis of large built-up structures. A general purpose 
nonlinear dynamic analysis module, DYCAST, is in final stages of 
development. These modules are not included in this document. 

The modules that analyze material nonlinearities alone employ 
the "initial strain" concept within an incremental procedure to 
account for the effect of plasticity and include the capability 
for cyclic plastic analysis. The "initial strain" approach does 
not require that the stiffness matrix be updated at each step in 
the analysis but rather, the effect of plasticity enters into the 
analysis as an effective load vector. The cyclic plastic' behavior 
is accounted for by implementing the Prager-Ziegler kinematic 
hardening theory. Geometric nonlinearities are included in the 
PLANS system programs by making use of an "updated" or convected 
coordinate approach, which requires the reformation of the stiff- 
ness matrix due to changes in the geometry and stress field during 
the incremental procedure. 

The theoretical foundations upon which the PLANS system pro- 
grams are based can be found in a companion volume, "PLANS - A Finite 
Element Program for Nonlinear Analysis of Structures, Volume I - 
Theoretical Manual," NASA Contractor Report NASA CR-2568. 

This present volume describes the input data preparation for 
the four principal modules of PLANS, BEND, HEX, REVBY, and OUT-OF- 
PLANE MG as well as the SATELLITE program for input checking. In 
keeping with the philosophy of the PLANS system, the description 
of the input for. each of the modules is in a self-contained sec- 
tion. Also presented is a section that shows the input data decks 
for a representative number of sample problems. 
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The input for all the modules begins with a title card that 
allows for any 80 character title (specified in columns 1-80) . 

This title serves as a page heading for subsequent computer out- 
put, The input data following this card is divided into a number 
of functional groups, each describing a specific type of input 
information. These groups are briefly described below and sche- 
matically shown in Fig. 1. Each input group must be read in the 
specific order shown in the figure. In general, each group is 
delineated with a specific section end card. This is the alpha- 
numeric SEND, left justified on the input card, in columns 1 through 
4 or a blank card. The input groups are as follows: 

I Title Card 

II Program Control Parameters and Options 

III Node Specification 

This section defines an allowable set of node 
point identification numbers 

IV Element Connectivity 

-Defines each element by specifying its type 
(i.e,, beam, triangle, ..., etc .), identifica- 
tion number and connecting node points 

V - VII Node Point Coordinates 

Defines the location of each node point in a 
global cartesian coordinate system 

VIII Node Point Single and Multipoint Constraints 
Defines boundary conditions and nodal con- 
straint equations 

IX Element Material and Section Properties 
Section properties include element thickness, 
area and moment of inertia where applicable. 


1-4 



ShND- 


APPLIED LOAD DISTRIBUTION 
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U.KMI'NT MATER I Al. AND SECTION PROPER'! I TS 
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NODE SPEC1 1' ICA I ION 


PROGRAL1 CONTROL PARAME fKRS 


PROBLEM TITLE CARD 


Fig. 1 Order of Groups of Input Data 


1-5 



Material properties include such quantities 
as Young’s modulus', Poisson's ratio, and 
quantities defining the plastic response such 
as yield stress and hardening parameters, 

X Applied Load 

Defines spatial distribution of applied load 

The last card in the deck is an alphanumeric STOP or EM) left 
justified on an input card in columns 1 through 4. STOP indicates 
that the job is complete and END indicates another problem deck 
follows . 

Some general rules have been used in designing the input. 

These rules are listed below: 

® Some of the data sections make use of a "keyword" 
of up to four characters, left justified in their 
appropriate field to specify an item of data. For 
example, cards specifying element connectivity for 
membrane triangles begin with TRIM, and material 
properties by MATl. All sections that make use of 
keywords end with a section end keyword SEND 

e Sections that specify data without keywords end with 
a blank card section delimeter 

* The data deck ends with one of two keywords. If END 
is used, the program reads a new problem data deck. 
If STOP is used, as will probably be the case for 
most problems, the job ends 

9 Two formats for input are used, E15.7 for floating 
point input (fields of 15) and 15 for integer or 
fixed point input (fields of 5) . The fixed point 
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(integer) data must be right justified. The 

* 

floating point data can be written in several 
forms. For example, 10.0 can be input as: 

10.0 any place in the field, or 

1.0E+01, 1.0E+1, 1.0E1, where the 

field is right justified 

There are a number of places in the program where 
applicable node points or elements must be speci- 
fied with a set of data. In these cases the nodes 
or members are specified by entering the appro- 
priate number on the input cards in fields of five. 
However, for this purpose the user can also utilize 
a shorthand form of the input. That is, specifying 
m and -n consecutively is the equivalent of the 
specification of nodes or elements m , nH-1 , ro-t-2 , 

. , . , n and specifying m, -p, and -n consecu- 
tively is the equivalent of the specification of 
nodes (elements) m, m+p, m+2p, mlkp where 

tn+kp is the highest integer of the form less than 
or equal to n. For example, the specification of 
nodes 1 through 100 is written as 1-100 and 
nodes 1, 3 , 5 , . . . , 99 as 1, -2, -99. This card 
input appears in fields of 5 (15 format) with 36 
items per card. Any number of continuation cards 
may be used. A blank 15 field ends the specification. 
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SECTION 2 


INSTRUCTIONS FOR USE OF BEND 

A Program for the Elastic-Plastic Analysis 
of Built-up Structures 



INSTRUCTIONS FOR USE OF BEND 


BEND is a finite element program to treat the elastic, elastic 
plastic or elastic-cyclic plastic response of arbitrary built-up 
thin walled structures where bending and membrane effects are 
equally important. The finite element library consists of the 
following elements: 

® Three -node uniform strain triangle 
8 Six-node linearly varying strain triangle 
8 Four- and five-node hybrid triangles to be used 
as transition elements between three- and six-node 
elements (see Fig. la) 

® Two-node uniform strain stringer 
® Three-node linearly varying strain stringer 

• Beam with various cross sections subjected to 
bending about two planes as well as torsion 

• Higher order triangular plate element with 
bending and membrane capability 

The program is capable of treating the elastic and elastic- 
ideally plastic, linear strain hardening, and nonlinear strain 
hardening behavior of orthotropic materials. Further, the kine- 
matic hardening theory of plasticity is used (Refs. 1-3) so that 
provision for cyclic loading conditions involving reversed plastic 
deformation is included. 

The input to the program is categorized in the following sec- 
tions : 

I. Problem Title FORMAT (20 A4) 

Any 80-character title describing the problem. 

H . NPNTC, NPRNT, IRESRT, NUTAP, INPRT FORMAT (515) 

0 < NPNTC < 63: 
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NPNTC is the sum of the following integers corresponding to 
the option desired. 

If NPNTC = 0 No intermediate printout 

= 1 Print the load vector 

= 2 Print element stiffness matrix 

= 16 Print each element stiffness matrix entry 
to be stacked with its stacking index 
= 32 Print the total stiffness matrix 

For example, if it is desired to print the load vector and 
the total stiffness matrix, NPNTC = 1 + 32 = 33. 

NPRNT: 

If <0, perform elastic analysis only 

If > 0, perform plastic analysis, printing output every 
NPRNT increments of load 

IRESR.T: (See Section XVIII. ) 

0, Elastic (and/or)- plastic run. Do not generate 
a restart tape 

1, Elastic (and/or) plastic (cyclic) run. Build a 
new restart tape 

2, Plastic run with elastic values from previously 
created restart tape 

3, Plastic restart run starting at some specified 
load level 

4, Cyclic restart run 

NUTAP : 

Applicable only if restarting from a restart tape 

0, No new tape written 

1, Complete new -tape created and additional 
restart data written 
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INPRT : 


0, Write restart tape only at P = PMAX (i.e., at 
maximum load) 

N, Write restart tape every N increments of load 

III . Node Specification (1615) 

This section defines an allowable set of external node point 
numbers. The maximum node number that can be used is 9999. The 
program uses this information in two ways. First to set up a table 
of allowable node points that is used to check all subsequent node 
point input. Secondly, the program converts each external node num- 
ber to an internal number consecutively in the order that the node 
appears on the input card. Consequently the order of the input of 
external node numbers is completely arbitrary and need not be in- 
creasing monotonically. In practice the node' numbers should be num- 
bered so as to minimize the bandwidth. Once the input is read the 
program operates with the internal numbers which are now numbered 
frpm 1 through the number of nodes in the model. In this manner 
tf|p node ordering and therefore the bandwidth of the stiffness 
maprix can be easily^ changed and nodes can be inserted or deleted 
by phanging the external node specification. 

The input is specified by entering the appropriate number on 
the input cards in fields of five. However, for this purpose the 
user can also utilize a shorthand form of the input. That is, spe- 
cifying m and -n consecutively is the equivalent of the speci- 
fication of nodes m, m + 1 , m + 2, ..., n and specifying m, -p, 
and -n consecutively is the equivalent of the spe^ifict i.on of 
nodes m, m + p, m + 2p, . . . , -m + kp where m + kp is the highest 
integer of the form less than or equal to n. For example, the spe- 
cification of nodes 1 through 100 is written as 1 - 100 and 
1, 3, 5, ..., 99 as 1-2-99. This card input appears in fields of 
5 (15 Format) with 16 items per card. Any number of continuation 

cards may be used. A zero or blank 15 field ends the specification. 
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TRIM - Triangular membrane element (Ref. 3) 

BEAM - Beam element 

STRG - Stringer element (Ref. 4) 

TRIP - Triangular plate bending element (with 
membrane effects) (Ref. 5) 


The first integer field designates the member number. The 
next integer fields designate the connecting nodes as follows: 

TRIM - The nodes for the triangular family of elements 
are specified around the perimeter beginning with a major (vertex) 
node, and followed by a minor (midside) node and then alternatively 
major and minor as shown in Fig. 1. The absence of a minor node must 
be indicated by a zero or blank field in the proper position. 

BEAM - Three node specifications are necessary for the 
beam element . 


Nodes i and j (Fig. 2) which designate the element end 
points and a third node k, defining the normal to the beam axis 
about which the section properties are specified. This additional 
node may be a node of the structural idealization or it may be a 
"fictitious node" specified just for the purpose of defining the 
beam section properties. This is shown in Fig. 2. 

STRG - Three node specifications are necessary for the 
stringer element. Nodes i and j , connecting the end points and if 
desired an additional node designating a midpoint node. This is shown 
in Fig. 3. A zero or blank for the midpoint node specifies a two-node 
stringer. The midside node is the third node specified. 
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TRIP 


Three nodes are used to specify the plate 
element i, j, k as shown in Fig. 4. 

If a planar structure is being analyzed using triangular plate 
and membrane elements, stress and strain output can be calculated - 
with respect to the global coordinate system rather than the element 
local system by specifying one of the following cards: 


GLXY - 

Structure 

in 

. the 

x-y 

plane 

GLXZ - 

Structure 

in 

the 

x-z 

plane 

GLYZ - 

Struc ture 

in 

the 

y-z 

plane 


SEND - Ends the section. 

V. X-Coordinates of Nodes FORMAT (E15.7. 1315) 

The x-coordinates of the nodes appearing in the 15 fields 
are set to the value in the E15.7 field. Any number of continua- 
tion cards may be used; their first fifteen columns are ignored. A 

zero or blank 15 field terminates the card scan for a given 
x-coordinate. A zer o or blank first 15 field (columns 16-20) on 
a noncontinuation card ends the section. Both shorthand repre- 
sentations of Section III are allowed. 

VI. Y-Coordinates of the Nodes. Same as Section V. 

VII. Z-Coordinates of the Nodes. Same as Section V. 



VIII. Boundary Conditions FORMAT ri2ll. 3X . 13TS. /. I5X. 1315 


The first twelve fields give the boundary conditions specifica- 
tions in the order: u, v, w, 0 , 0 , 0 , e . e , e , k , k , k ; 

* ’ * x’ y z’ x y xy x y xy’ 

where u, v, w are global displacements in the x, y, z direc- 
tions, respectively; 6^, 0 , 0^ are global rotations with re- 
spect to the x, y, z axes; e , e are membrane components 

of the local membrane strain tensor; k , k , k are components 

x y xy r 

of the local curvature tensor. The last six fields for the local 
membrane strains and curvatures (higher order degrees of freedom) 
are only used if a plate element is attached to the node. Both ten- 
sors are calculated in the local coordinate system (see Fig. 4) of 
the plate element. This local coordinate system is defined by the 
first plate element connecting that node that appears in the - input 
stream. For example, 


TRIP 

3 

6 

3 

•2 " 


Node 3 is not listed pre- 

TRIP 

4 

3 

6 

7 

)■ 

viously in the element 
input 

TRIP 

5 

7 

4 

3 







J 




The tensor quantities at node 3 are all defined in the local coor- 
dinate system of element 3. See Appendix IV for examples of plate 
boundary conditions. 


The rotational degrees of freedom 9 6 , 6 ^ are applicable 

only for beam and plate elements. If nonzero values are given to 
rotational or higher order degrees of freedom at a node and beam 
and/or plate elements are not attached to the node then a singular 
stiffness matrix will result. 
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Zero denotes a fixed degree of freedom 

One denotes a free degree of freedom 

Two will result in the application of a unit gen- 
eralized displacement-, or a corresponding card 
may be included in Section IX (dependent degrees 
of freedom) designating the magnitude of the 
applied generalized displacement. 

The 1315 fields give the applicable nodes for the boundary 
condition specifications, with both shorthand notations of Sec- 
tion III permitted. Any number of continuation cards may be used 
for a given specification. However, only the 1315 fields are 
used on a continuation card. A zero or blank 15 field- terminates 
the card scan for a given boundary condition specification . Note: 

If the last field of a card (columns 76-80) is the last specifica- 
tion, an additional blank card (continuation card) must follow. A 
zero or blank first 15 field (columns 16-20) on a noncontinuation 
card ends the section. If a node's boundary conditions are not speci- 
fied in this section, all the degrees of freedom are assumed to be free. 
To change this default condition, the first card of this section should 
be set to the desired default (i.e., in the absence of beams and plates, 
111000000000) with all nodes used in the problem specified. In the ab- 
sence of plate elements the default is automatically 111111000000. 

IX. Dependent Degrees of Freedom FORMAT (215, 2(5X, 215, E15.7 )) 

This section designates the input for both single and multipoint 
constraints as well as applied displacements of the form: 

1) 6. = ol 6 . - +a 0 6. 0 4- •••a 5. 

i T. jl 2 j2 n jn 

where 5^ is a dependent degree of freedom, 

6.. **• F>. are independent degrees of freedom, 
and a. , a„ • • • a are coefficients . 

L L T1 
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2 ) 


Rotation of displacements at a node 


5. 

l 


cr 6. + a^5. + a 0 5, 
1 i 2 j 3 k 


o. = P 1 6. + & 2 S. + P 3 6 
6 k - ■y 1 B i + y 2 l } + r 3 \ 


where the a, (3, ■y’s are the direction cosines of 
the rotation, 6., 5.., ft are the displacements with 

1 /v 

respect to the original global directions and 6., 

5., 5, are the components of displacements at the 
J 

node with respect to the new coordinate axes . An 
example of this capability is given in Appendix I. 


3) Applied generali 2 ed displacement 


5 . = a 
x 


where the coefficient a is the applied generalized 
displacement . 


The first two-fields designate a node number and a degree of 
freedom (i.e., 1-12). The dependency is defined in the following 
three fields. The two integer fields designate the node number 
and degree of freedom number and the coefficient by the floating 
point field. If there is another dependency for the node, it is 
designated in a similar fashion in the next three fields. Any 
number of continuation cards can be used with the first too fields 
blank. The section is ended by a blank or zero in the third in- 
teger. field (blank card) . Examples of the use of multipoint con- 
straints are in Appendix I. 
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X. 


Material and. Section Properties 


The format for this input is dependent upon the member type. 

• Each type of input begins with a word of up to four characters. 

MATl Material properties for plane stress, necessary 
with membrane triangles (TRIM) and bending 
triangles (TRIP) . For use with initially iso- 
tropic materials which exhibit perfectly plastic 
or strain hardening behavior and initially 
orthotropic materials that exhibit perfectly 
plastic behavior. 

MAT2 Material properties for plane stress, necessary 
with membrane triangles (TRIM) and bending tri- 
angles (TRIP) with initially isotropic or ortho- 
tropic behavior that require orthotropic kinematic 
hardening theory 

MBM Material properties for a beam element 

MSTG Material and section properties for a 
stringer 

THIK Member thickness for triangular membrane (TRIM) 
and plate (TRIP) elements 

MBET Angle between local axes and principal direc- 
tions of orthotropy for TRIM and TRIP elements 
(see Figs, lb and 4) 

NLRfcJ Layer and layer output information for TRIP 
elements 

SREG Beam section properties for a solid rectangular 
section ’ 
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SCIR 

Beam section properties 
section 

for 

ZSEC 

Beam section properties 

for 

ISEC 

Beam section properties 

for 

HCIR 

Beam section properties 
section 

for 

HREC 

Beam section properties 
section 

for 

LSEC 

Beam section properties 

for 

TSEC 

Beam section properties 

for 

CSEC 

SEND 

Beam section properties 
Ends the section 

for 


a solid circular 

a Z-section 
an I- section 

a hollow circular 

a hollow rectangular 

an L-section 
a T-section 
a Channel section 
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MA.T1 - Plane Stress Material Properties - FORMAT (A4, IX, 5E15.7, /, 

2E15.7, /, 4E15.7, /, 
5E15.7, (, (1615)) 

The first four cards specify material properties, as fol- 
lows : 

CARD 1: MATl 

EONE = Young’s modulus in principal property axis (1) 
ETWO = Young's modulus in principal property axis (2) 
BETA - No longer used. Set equal to zero. 

GONTO =* Shear modulus in (1) - (2) principal property 
plane 

VONTO - Poisson's ratio, V 12 


Note: 




^2 E x °1 + E 2 


7 


12 



CARD 2: TALF-1 - Coefficient of thermal expansion in 

1- axis direction, 

TALF-2. = Coefficient of thermal expansion in 

2- axis direction, 

CARD 3: SIGOX - Yield stress in principal 1-direction 

SIGOY = Yield stress in principal 2-direction 

SIGOZ « Yield stress in principal 3-direction 

SIGXY = Shear yield stress in principal 1-2 plane 
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CARD 4: RMOSN = If RMOSS ^ 0; RMOSN = n, the shape 

parameter used in Ramberg-Osgood repre- 
sentation of stress -strain behavior 

If RMOSS = 0; RMOSN = a the slope of 
the linear strain-hardening stress-strain 
representation, i.e., a = E^/E where 
is the tangent modulus 

RMOSS = If RMOSS ^ 0; RMOSS = Ramberg-Osgood 
parameter ^ (see Ref. 7) 

Note 1: If RMOSN = 0 and RMOSS = 0, 

the material for the element (s) 
is assumed to be elastic-ideally 
plastic 


RMOSE = Ramberg-Osgood parameter E (Young's modulus) 

YLDST - Yield stress ,in tension 

YLDSC = Yield stress in compression 

Input for the yield-stress in tension and compression has been main- 
tained in order to accommodate materials that exhibit initial aniso- 
tropic plastic behavior. In this case an initial translation of the 
yield surface is made consistent with the kinematic hardening theory. 


Note: Only initially isotropic materials can be treated when con- 

sidering linear or nonlinear strain hardening using this 
property card. 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be used 
for a given specification. A zero or blank 15 field ends each member 
listing . 


Isakson, G. , Armen, H. Jr. , and Pifko, A. , "Discrete-Element Methods 
for the Plastic Analysis of Structures," NASA. Contractor's Report 
NASA CR-803, October 1967. 
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MAT2 - Plane Stress Material Properties - FORMAT (A4, IX, 5E15.7, /. 

2E15.7, /, 4E15.7, /, 
5E15.7, I, E15.7, /, 
(1615)) 

The first five cards specify material properties, as fol- 
lows : 


CARD 1: MAT 2 

EONE - Young's modulus in principal property axis (1) 
ETWO = Young ' s modulus in principal property axis (2) 
BETA = No longer used. Set equal to zero. 

GONTO = Shear modulus in (1) - (2) principal property 
plane 

VONTO = Poisson's ratio, V 12 

C 1 v 12 

Note: e l "5 ^'e^~ 9 2 





Y 12 G 12 T 12 


CARD 2: TALF-1 - Coefficient of thermal expansion in 1-axis 

direction, 

TALF-2 = Coefficient of thermal expansion in 2-axis 
direction, 

CARD 3: SIGOX = Yield stress in principal 1-direction 

SIGOY = Yield stress in principal 2 -direction 

SIGOZ = Yield stress in principal 3-direction 

SIGXY = Shear yield stress in principal 1-2 plane 
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CARD 4: RMONX = If RMOSX $ 0-; RMONX = n ; the shape 

parameter used in the Ramberg u dsgood 
representation. of the com- 

ponent of the stresses train Behavior. 


If RMOSX = 6; RMONX = a the slope of 
the linear strain -hardening stress -strain 


repf eseiif at ion iri the 

— E t1 

“i = IT' whe ” e ti 


l-direction, i;e 
is the tangent 


9 


modulus of the 


T 1 


curve . 


RMOSX 

If RMOSX 4 0; 

RMOSX = Ramberg-Osgood 


parameter o ^ ^ 

in principal 

1-direction 

RMONY 

Same as 

RMONX 

in principal 

2-direction 

RMOSY 

Same as' 

RMOSX 

in principal 

2-direction 

RMONXY 

Same as 

RMONX 

for shear component in 


1-2 plane; “ 7^ stress-strain 

curve 


CARD St 


RMOSX-Y 


Safiie as RMOSX 
i-2 plane; 


for shear component in 
- 7-^2 stress-strain 


curve 


Note: No component of the stress -strain 

curves may Be ideally plastic > 
i . e . , RMONX or RMONY or 
RMONXY = 0 is not allowed.- If 
either is zero all must Be zero 
and a MATl card should be used. 

The Ramberg-Osgood representation of 
the stress -strain curve is given by 
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a 

- — + 


3 a 
7E 


(“) 


n-1 


0.7 


7 = 




n-X 


Succeeding cards give applicable members; both shorthands of 
Section III are permitted. Any number of continuation cards may 
be used for a given specification. A zero or blank 15 field 
ends each member listing. 
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MM - Beam Material Properties - FORMAT (A4, IX, 5E15.7, /, E15.7, / , 

(1615)) 


CARD 1: MBM 

E = Young 1 s modulus 

ANU = Poisson’s ratio 

RMOSS = If not equal to zero, RMOSS equals 
Ramberg- Osgood parameter, cTq -j 

RMOSN = If RMOSS 0; RMOSN = n, the shape parame- 
ter used in Ratriberg-Osgood representation of 
stress -strain behavior. If RMOSS = 0; 

RMOSN = a, the slope of the linear strain- 
hardening stress-strain representation, i.e., 

a =5 E m /E where E„ is the tangent modulus 

r i 

YLDST = Yield stress 

CARD 2: TALF = Coefficient of thermal expansion 


Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing. 
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MSTG - Stringer Properties - FORMAT (A4, IX, 5E15.7, / E15.7 /, (1615)) 

CARD 1: .MSTG 

E = Young's modulus 

A = Cross-sectional area 

RMOSN = If RMOSS ± 0; RMOSN - n, the shape 

parameter used in Ramberg -Osgood repre- 
sentation of stress-strain behavior. 

If RMOSS = 0; RMOSN = a, the slope of 
the linear strain -hardening stress -strain 
representation, i.e., a - E^/E where E^ 
is the tangent modulus 

RMOSS = If RMOSS $ 0; RMOSS « Ramberg -Osgood 
parameter Oq y 

Note 1: If RMOSN =. 0 and RMOSS = 0, 

the material for the element (s) 
is assumed to be elastic-ideally 

t 

plastic 

YLDST = Yield stress 

CARD 2: TALF = Coefficient of thermal expansion 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing . 
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THIK - Element Thicknesses - FORMAT (A4, IX, E15.7, /, (1615)) 

(Necessary with membrane elements 
TRIM) and plate elements (TRIP)) 

CARD 1: THIK 

THICK *= Element thickness 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing . 


MBET - Orientation of Axes of Material Anisotropy - FORMAT (A4, IX, 

E15.7, I > (1615)) 


CARD 1: MBET 

BETF - Angle (3 in degrees between local x-axis 
and principal 1-axis for material ortho- 
tropy. See Figs, lb and 4. Only appli- 
cable for TRIM and TRIP elements 


Note: This card is an optional card. The 

default is BETF = 0.0 for all members. 
It should be used with initially ortho- 
tropic materials or those exhibiting 
ortho tropic strain -hardening character- 
istics . 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing.. 
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NLRS - Layer and Print Option for TRIP Elements - FORMAT (A4, IX, /, 

215) 

CARD 1 : NLRS' 

CARD 2: NLRS = Number of layers through the thickness at 

which stresses are calculated for TRIP 
elements . NLRS must be an even number 
•< 20 

NPLT = Print out option: 

= 0; Print stresses, strains, etc. only at 
top and bottom surfaces 
- 1; Print stresses, strains, etc. at each 
layer through the thickness as well. 

Note: This is an optional card. The 

defaults are NLRS = 10 and 
NPLT = 0 . If one is changed 
both should be specified. This 
card controls option for all 
TRIP elements . 
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Beam Section Properties - FORMAT (A4, IX, 5E15.7, l , 3E15.7, /, 

formats for Cards 3 and 4 (see 
below), /, (1615)) 

The cards specifying beam section properties all start with 
the following information: 


CARD 1: 


SRECh 
SCIR 
1SEC 
ZSEC 


Defined at the beginning of the section. 
HCIR \ The different cross sections are shown 


HREC 
LSEC 
TSEC 
CSEC j 
A 

fyy 

Izz 
Iyz 
J ' 


in Tables 1 and 2. 


= cross-sectional area 
= moment of inertia about 
= moment of inertia about 
= product of inertia 
= torsional rigidity 


y axis (see Fig. 
z axis (see Fig. 


5) 

5) 


If the quantities, A, I , ^ zz > ^ y Z are input as zero or 
blank then these quantities are calculated automatically. 



Fig. 5 Definition of Coordinate Axes in 
Cross Section of Beam Element 
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CARD 2 : YO 


ZO 

P 


- Eccentricity of attachment point in the y Q 
direction 

= Eccentricity of attachment point in the z q 
direction 

= angle defining the transformation of the 
Iy, Iz, Iyz to another reference axis 
(see Fig. 5) 


Additional cards are required according to which section is speci- 
fied. The notation for each section is shown in Table 1. 


SREC - FORMAT (2E15.7) 

CARD 3: A = Width 

B = Depth 

SCIR - FORMAT (E15.7) 

CARD 3 : • R = Radius 

ZSEC, ISEC, and CSEC - FORMAT (3E15.7, /, 3E15.7) 

CARD 3 : A1 = Dimension of upper flange (parallel to y-axis) 

A2 = Dimension of web (parallel to z-axis) 

A3 = Dimension of lower flange (parallel to y-axis) 


CARD 4: Tl = Thickness of upper flange 

T2 - Thickness of web 

T3 = Thickness of lower flange 

HCIR - FORMAT (2E15.7) 

CARD 3: R = Outer radius 

T = Thickness 
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HREC, LSEC, and TSEC - FORMAT (4E15.7) 


CARD 3 : Al 

A2 
T1 
T2 


= Width (parallel to y-axis) 

- Depth (parallel to z-axis) 

= Thickness of upper and/or lower flanges 
= Thickness of vertical webs 


Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing, 

SEND Ends the section. 
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Section 

Keyword 

Input 

Solid 

Rectangle 

SREC 

! 

i 

i 

i 

Area, moments of inertia, 1,1 , I ' , 

’ ’ yy’ zz’ yz’ 

J eccentricities, y , z ; transformation 
' J o o 

angle Section width and height, a, b. 

fSolid 

Circular 

SCIR 

Area, moments of inertia, I ,1,1, 

' ’ yy’ zz- yz' 

J eccentricities, y , z , transformation 
angle p . Section radius, r. 

Hollow 

Circular 

Section 

HCIR , 

. 

Area, moment of inertia, I , I , I , 

yy zz yz 

J eccentricities, y , z , transformation 
_ o o 

angle f3. Outer radius, r, thickness, t. 
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TABLE 1 SECTION INPUT -AND GEOMETRY (CONT) 


Section 

Keyword 

Input 

Hollow 
Rectangular ■ 
Section 

HREC 

Area, moment of inertia, I , I , I , 

’ 1 yy’ zz yz ■ 

J eccentricities, y , z , transformation 
angle 3. Width, >a^, depth, thickness of 

upper and lowerj flanges, t^, thickness of 
vertical webs, t^ 

L- Sect ion 

LSEC 

Area, moment of inertia, I, I , I , 

yy y z 

J eccentricities, y , z , transformation 
angle 3. Dimension of flange and web apa£ 
thickness flange, and web, t^. 

T- Sect ion 

TSEC 

Same as for L- Section 



iffm 

. — A i — 
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TABLE 1 SECTION INPUT AND GEOMETRY (CONT) 


Section 

Keyword 

Input 

...... .... - 



Z-Section 

ZSEC 

Area, moments of inertia, I, I , I , 

y y zz 

J eccentricities, y , z , transformation 

o o 

angle p. Dimension of upper and lower 
flange, a^, a^j Dimension of web, a 2 » 
thickness of upper and lower flange, t^, t^, 
thickness of web, t 2 
NOTE: a^ = a^; = t^ 


I 2 

1 T 2 Y t 

"tP 

a 2 

j 



~ 

I-Section 

I SEC 

Same as Z-section 
a l ^ a 3 ; t l ^ t 3 

K A 1 


I T * 

i 


H" T 3 * ’ 

K- 

* 

A 3' — H 

Channel 

CSEC 

Same as Z-Section 

t 

7 

Section 



r 




a l = a 3 * C 1 = t 3 

TF 

± S. tJ 

T i 





►y 




I • 

^2 

T 3 



■ - - 

H“ A 3->-| j 




2Z- 


TABLE 2 NUMBER AND LOCATION OF STRE 


Keyword | Number 
Stress 
Points 


Location of 
Stress Points 


Solid 

Circular 


± y 

0.430568a 
0.430568a 
0. 169991a 
o! 169991a 


+ z 

0.430568b 
0. 169991b 
0.430568b 
0.169991b 



Location of 
Stress Points 














6Z- 


N> 


Section 


Hollow 

Circular 

Section 


Hollow 

Rectangular 

Section 


Keyword 


HCIR 


HREC 


Number 

Stress 

Points 


8 


12 


1 

2 

3 

4 

5 

6 

7 

8 


1 

3 

7 


9 


lND location of STRESS POINTS (CONT) 


Location of 
Stress Points 


Location of 
Stress Points 


0. 77459a, 


0. 5a. 


0.5a, 


352.8 


323.4 

36.6 

274.6 

85.4 

213.0 

147.0 


0. 5a,. 


0. 5a r 


0. 77459a, 



Shear Center at Centroid 



a l “ a l " C 2 


Shear Center at Centroid 


a 2 a 2 ” t l 
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TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 


Section 


Keyword Number 
Stress 
Points 


|L~ Sect ion 


Location of 
Stress Points 


5 | 

2 3 1 

-* ft 2 


Axes at Shear Center 



Axes at Shear Center 



TC- 


TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 




Z£- 


TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 


Section Keyword Number 

Stress 



Location of 
Stress Points 



y 

z 

1 

0.88730a, +e 
1 s 

0.5a, 

2 

0. 11271a, + e 
1 s 

0.5 a 

3 

0.5000a, + e 
1 s 

0. 5a 

4 

0. 88730a, + e 
1 s 

-0.5a 

5 

0. 112 71a, + e 
1 s 

~0.5a 

6 

0.5000a, + e 
1 s 

-0.5a 


0. 77459a. 


-0.77459a. 


e - distance from shear center to 

g 

web center line 


e = 
s 


1/2 (a L - t 2 ) 

. - , J t 2 (a 2 “ t l 


1 + 1/6 


t 9 (a, - t,)\ 1 a 2 a 2 _t l 


t i (a r " t 2 ) i 


NOTE: t x = t 3 


a i “ a 3 


Location of 
Stress Points 



Axes at Shear Center 
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TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 


ro 









XI.. Applied Loads 

Many different types of loading are currently admissible — 
concentrated loads at nodes, distributed loads on an edge of a 
triangular member, distributed surface loads, distributed line 
loads on beams, and thermal loads. Each loading situation is 
designated by one of the following four-character key words: 

CONC (Concentrated force or moment at a node) 

TRIA (Distributed edge load, TRIM elements) 

SURF (Distributed surface loads, TRIP element) 

EDG1 (Distributed edge load, side 1, TRIP element) 

EDG2 (Distributed edge load, side 2, TRIP element) 

EDG3 (Distributed edge load, side 3, TRIP element) 

BMLO (Distributed line load, beam element) 

TMPU (Temperatures at upper surface for TRIP 

element or temperatures for TRIM element) 

TMPM (Temperatures at middle surface for TRIP 
element) 

TMPL (Temperatures at lower surface for TRIP 
element) 

SEND (Section end) 

Note: The keywords for the section on applied loads 

are used as section headings. The keyword 
appearing on an input card designates that the 
input to follow is associated with a particular 
type of applied loading. A blank card (as de- 
scribed in succeeding sections) delimits the 
input and indicates that the next card contains 
a different keyword. For example, CONC indi- 
cates input that follows is for concentrated 
loads and TRIA; distributed loads on an edge of 
a membrane triangle. Thus the input is as follows: 
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CONC 


data 

for 

concentrated 

loads 

blank card 

TRIA 

data 

for 

distributed 

loads 

blank card 

SEND 

It should be noted that the keywords can appear in 
any order and may be specified more than one time 
in the applied load section. 

TRIA - Distributed edge load (TRIM elements) - FORMAT (A4, /, 415, 

3E15.7, /, 15X, 

3E15.7) 

Each card gives the load components applied at a member side 
(in units of force) as follows: 

CARD 1: TRIA 

CARD 2: First 15 field: Number of node (m) 

Second 15 field: Number of minor node on the side, 

aero, or blank if there is no 
minor node 
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.Number of other end point node (n) 
NODR — Reference node number 
x load component at node m 
y load component -at node m 
z load component at node m 

x load component at node n 
y load component at node n 
z load component at node n 

Note: If NODR = 0 then x, y, 

z are global directions . 

If NODR ± 0 then the 
x-direction is parallel 
to the element edge speci- 
fied by m - n; the 
y-direction is perpendicu- 
lar to x in the plane 
defined by the three points, 
m s n, NODR; the z-axis is 
perpendicular to the x-y 
plane. 

The program allows for a linearly varying distributed load from node 
to^node. As many load components as desired may be specified. A blank 
card (i,«e.., zero or blank first 15 field) ends the specification of 
distributed edge loads for TRIM -elements . 


Third 15 field: 
Fourth .15 field-: 
First E15.. 7 field: 
Second E15.7 field: 
Third E15.7 

CARD 3: First E-15-.7 field: 

Second E15.7 field: 
Third E15.7 field: 
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CONC - Concentrated loads - FORMAT - (A4, /, 15, 3E15.7, /, 5X, 

3E15.7) 

Each card gives the load components in global directions at a 
specified node, 

CARD 1: CONC 

CARD 2: 15 field: . Node number 

First E15.7 field: Force component F 

A 

Second E15.7 field: Force component F^ 

Third E15.7 field: Force component F> 

A 

CARD 3: First E15.7 field Moment component M 

Second E15.7 field Moment component M 
Third E15.7 field Moment component M 

A 

A blank card (i.e., zero or blank first 15 field) ends the spe- 
cification of concentrated loads - 

SURF - Distributed Surface Loads (TRIP element) - FORMAT (A4, IX, 

3E15.7, /, 3E15.7, /, 
3E15.7, /, (1615)) 

CARD 1: SURF 

PXX Value of distributed surface load in force/ 
unit area in local element x-direction at 
node i of element (see Fig. 4). 

PYI Value of distributed surface load in force/ 
unit surface area in local y-direction at 
node i of element 

PZI Value of distributed surface load in force/ 
unit surface area in local z -direction at 
node i of element 
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CARD 2: 


CARD 3: 


PXJ Same as PXI but at node j 
PYJ Same as PYI but at node j 
PZJ Same as PZI but at node j 

PXK Same as PXI but at node k 
PYK Same as PYI but at node k 
PZK Same as PZI but at node k 

Note: A linear variation of the surface loads 
in the plane of the element is assumed. 


Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing . 


BMLO - Distributed Line Load, Beam Element - FORMAT (A4, /, 4E15.7, /, 

(1615)) 


Card 1: 
CARD 2; 


BMLO 

PYI 

Force /unit length in 

local 

y-direction 

(see 

PYJ 

Fig. 5) at node i 
Force/unit length in 

local 

y-direction 

at 

PZI 

node j 

Force /unit length in 

local 

z -direction 

at 

PZJ 

node i 

Force/unit length in 

local 

z -direction 

at 


node j 

Note: A linear variation of 

the distributed 


load between nodes is assumed. 


Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing. A blank card ends this section. 
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EDGl - Distributed Edge Load, Side 1, TRIP Element - FORMAT (A4, IX, 

3E15.7, /, 3E15.7, 
/, (1615)) 

This card inputs the distributed edge or line loads along side 1 
of a TRIP element. Side 1 is the side connecting nodes i and j 
(see Fig . 4) . 

CARD 1: EDG1 

PXI Local x-component of force /unit length at 
i^ node 

PYI Local y-component of force/unit length at 
i**k node 

PZI Local z -component of force /unit length at 
i fc k node 


CARD 2: PXJ Same as PXI but at node j 

PYJ Same as PYI but at node j 

PZJ Same as PZI but at node j 

Note: A linear variation of the edge loads 

is assumed between nodes i and j . 


Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 


used for a given specification. A zero or blank 15 field ends 
each member listing. 


EDG2 - Distributed Edge Load, Side 2, 


TRIP Element - FORMAT (A4, IX, 

3E15.7, /, 3E13.7 
/, (1615)) . 


This section inputs the distributed edge or line loads along 
side 2 of a TRIP element. Side 2 is the side connecting nodes j 
and k (see Fig. 4) . All input is the same as in the EDGl sec- 
tion with EDG2 replacing EDGl as the alphameric clue word and 
" j " replacing "i" and "k" replacing "j" in the description. 


2-39 



EDG3 - Dis.tr rbuted Edge .Load, Side 3.,. TRIP 'Element - FORMAT (A4, IX, 

3EI5.7, /, 3E15.7, 
l, (1615)) 

.This section inputs the distributed ;edge or line loads along 
side 3 of *a TRIP element. :Side 3 is the side connecting nodes k 
and i (see .’Fig. 4) . All input is the ;same ’as in the : EDG1 sec- 
tion with EDG3 replacing ED’Gl as the alphameric 'due word and "k" 
replacing "1" and "I" replacing "k" in the description. 
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TMPU - Nodal Temperatures At Upper Surface for TRIP Element or 

Nodal Temperatures for TRIM Element - FORMAT (A4, /, E15.7, 

1315, /, (15X, 1315)) 

This section is used to input temperatures at nodes of TRIM 
elements and at the upper surface (i.e., z - + h/2, where h is 
the plate thickness) of TRIP plate elements. For TRIP elements a 
parabolic temperature distribution through the thickness is assumed. 
To input temperatures at the lower surface and middle surface use 
the TMPL and TMPM cards, respectively. 

CARD 1: TMPU 

CARD 2: E15.7 field Temperature at node (TRIM) 

or temperature at upper surface 
of node (TRIP) 

15 Integer fields Applicable nodes 

The temperatures of the nodes appearing in the 15 fields are 
set to the value in the E15.7 field. Any number of continuation 
cards may be used; their first fifteen columns are ignored. A zero 
or blank 15 field terminates the card scan for a given temperature. 
A zero or blank first 15 field (columns 16 -20) on a noncontinua- 
tion card ends the TMPU input . Both shorthand representations of 
Section III are allowed. 


Note: For the TRIM element the 

nodal temperatures are 
averaged to obtain one 
element temperature . For 
the TRIP element a linear 
in-plane variation of the 
temperatures from node to 
node is assumed. 
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TMPM - Nodal Temperatures at Middle Surface for TRIP Element - 

FORMAT (A4, /, E15.7, 1315, /, 
(15X, 1315)) 

This section is used to input temperatures at nodes at the 
middle surface (z = 0) of TRIP plate elements. For TRIP elements 
a parabolic temperature distribution through the thickness is as- 
sumed. To input temperatures at the upper and lower surface use 
the TMPU and TMPL cards, respectively. 

The input formats and information for the card' are identical 
to the TMPU card with TMPM replacing TMPU. 

t 

TMPL - Nodal Temperatures at Lower Surface for TRIP Element - 

FORMAT (A4, /, E15.7, 1315, /, 
(15X, 1315)) 

This section is used to input temperatures at nodes at the 
lower surface (z = - h/2, h is the plate thickness) of TRIP 
elements. For TRIP elements a parabolic temperature distribution 
through the thickness is assumed. To input temperatures at the 
upper and middle surface use the TMPU and TMPM cards, respectively. 

The input formats and information for this card are identical 
to the TMPU card with TMPL replacing TMPU. 

A card with SEND in the first four columns ends the section 
for applied loading . 

XII. Members to be Printed - FORMAT (1615) 

Specifies the members whose strains and stresses are to be 
printed. Both shorthands of Section III are allowed. A blank card 
or card with only zero entries ends the section. 
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XIII. Nodes to Be Printed - FORMAT (1615) 


Nodes whose displacements are to be printed as per Sec- 
tion XII. 

XIV. Parameters for Plastic Analysis - FORMAT (2E15.7) 

PMAX - Maximum load to be applied for this half 
cycle 

PPCT = Load increment as a decimal multiple of 
yield load. 

XV. Parameters for Succeeding Load Cycles - FORMAT (I5~, 2E15 . 7) 

NPRNT = If equal to zero, no additional load cycle 
and end of problem. If greater than zero, 
print output every NPRNT increments. 

PMAX = Maximum load for this load cycle 

PPCT *= Load increment as percentage of yield load. 

Note: PPCT should be equal to one -half the 

original value since the plastic 
range is twice the elastic range for 
initial loading. 

XVI. Change Plastic Material Properties for New Load Cycle 

The input for this section begins with a four-charactep word 
as follows: 

MAT1 New material properties for plane stress 
membrane triangles (TRIM) or plate (TRIP) 
elements (see Section X for applicability) 

MAT 2 New material properties for orthotropic 

strain-hardening materials (see Section X for 
applicability 

MBM New material properties for beam elements 
MSTG New material properties for a stringer 
element 



MAT1 - FORMAT (A4, IX, 15, /, 4E15.7, } , 5E15.7) 


CARD 1: MAT1 

IMEM One member number from the group that is 
to be changed 


CARD 2: 


SIGOX s 
SIGOY 
SIGOZ { 


Defined in Section X 


SIGXY J 


CARD 3: 


RMOSN-\ 

RMOSS 


RMOSE > Defined in Section X 


YLDST 

YLDSC- 


One set of three cards is required for each MATl group to be 
changed. A MATl group may be changed to a MAT2 group after a 
half load cycle. 


MAT2 -FORMAT (A4, IX, 15, /, 4E15.7, /, 5E15.7, /, E15.7) 


CARD 1: MAT2 

IMEM One member number from the group that is to 
be changed 

CARD 2: SI GOX" 

SIGOY l 

) Defined m Section X 
SIGOZ f 

SIGXY J 


CARD 3: RMONX b 

RMOSX 


RMOMY ; Defined in Section X 

RMOSY 

RMONXY-' 
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CARD 4: RMOSXY Defined in Section X 

One set of four cards is required for each MAT2 group to be 
changed. A MAT2 group may be changed to a MAT1 group after a 
half load cycle. 

MBM or MSTG - FORMAT (A4, IX, 15, /, 3E15.7) 

CARD 1: MBM 

or 

MSTG 

IMEM (One member number from the group that is 
to be changed) 

CARD 2: RMOSN (Ramberg-Osgood parameters defined in Section X) 

RMOSS 
YLDST 

The section is ended with a section end card, SEND. 

Note: If there is no change in material 

properties, a SEND card must appear 
in the appropriate place in the 
input stream. 

XVII. Problem End - FORMAT (A4) 

Each problem's input must be ended with a card reading "ENDb" 
where "b n denotes a blank, in columns 1-4. The last problem- in 
a run should end with a card reading "STOP" in column 1-4 instead 
of "ENDb." 
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XVIII. Restarting a Problem 

The initial restart tape is created on Unit 2*1., Subsequent re- 
start jobs mount the restart tape- as- Unit 2-1 and' if desired .(NUTAP > 0) 
copies and continues the restart tape- on Unit 22. 

CARD; 1: Problem, Title - .FORMAT' (20A4.) 

As in. Section I 

CARD- 2:- NPNTG , NPRNT, IRESRT, NUTAP, INPRT ~ FORMAT (515), 

As in. Section II. Here NPNTC is ignored’. 

If IRESRT is equal to 2 

Restart from an elastic critical load using a- 
previously created restart tape 

CARD- 3’: PMAX', PPCT - FORMAT (2E15, 7) 

PMAX = Maximum, load to be applied- for this half cycle 

PPCT = .Load increment as percentage of yield’ load 

Input is then continued from- Section XV 

If IRESRT is eq,ua,l to 3 

Continuation from, some previously plastic load. The 
j.ob can be restarted from a. previous maximum load, or 
from any intermediate load level. 

CARD- 3: PMAX,, NRSRT - FORMAT (-E15.0:, 15) 

If starting from a previous maximum load' NRSRT can 
be ignored,. If starting from some- intermediate 
load level, NRSRT is a unique number obtained as 
output from the j-.ob- that generated the restart tape. 

This output is of the form-i 
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* * * * * 

PLASTIC ANALYSIS VALUES FOR RESTART TAPE HAVE BEEN 
WRITTEN AT P = l. 5000000E+03 FOR PRINT INTERVAL 30 

it it it it it 

In this case NRSRT is equal to 30 and the load is 
not necessarily the last maximum value. 

Note; For cases when the job ends abnormally (for example when 
time or lines are exceeded) the job can be restarted with 
a value for NRSRT obtained from the output. In this case 
the last value of NRSRT should not be used. This will 
insure that a complete set of data will be accessed from 
the restart tape. 

If IRESRT is equal to 4 

Restart from a previous maximum load level 
followed by a reversal of a previous load 
distribution (CYCLE). A reversal of a 
previous load distribution can be placed 
on the structure from any previously gener- 
ated and saved state of stress. 

CARD 3: PMAX, NRSRT - FORMAT (E15.0, 15) 

If starting from a previous maximum load NRSRT can 
he ignored. If smarting from some intermediate load 
Level, NRSRT is a unique number obtained as output 
from the job that generated the restart tape. Thi- 
output has the form shown in the preceding dis- 
cussion for IRESRT equal to three. The value of 
PMAX is ignored . 

Succeeding cards follow as in Sections XV through XVII. 
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XIX. Alternate Input Method 

All or part of the input may be read from a tape. This tape is 
mounted as Unit 23 with the appropriate job control cards. The tape 
is accessed by specifying the keyword, TAPE (FORMAT A4) , in the proper 
place in the card input stream. Input will then be read from the tape 
until a SEND (A4) is encountered on the tape. This returns the input 
reader to the card input stream. 

Note: All SEND formats used to end 

sections' must be in card form 
as the SEND encountered on the 
tape only sends the reader back 
to the card input stream. Also 
blank cards which end sections 
are best not put on the tape but 
rather on cards as this allows 
the user to add card input to 
tape input, i.e., if a blank 
ends a section and is encountered 
on the tape format then the sec- 
tion ends and no new data may be 
added, but if the blank does not 
occur on tape the section is 
still open and can be closed with 
a blank card in the card input 
stream. 

Example : 

Input for member connectivity. Section IV, and coordinates have 
been generated on a tape in the prescribed manner. 
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Input cards are as follows: 

TITLE 

0 5 

1 - 100 

TAPE (reads data from a tape) 

TRIM 

(additional elements if desired) 

* 

SEND 

TAPE (reads X-coordinates from tape) 

(Blank card to end section for X-coordinates) 
TAPE (reads Y-coordinates from tape) 

(Blank card to end section for Y-coordinates) 
5.0 (card input for Z -coordinates) 

(Blank card ends section for Z-coordinates) 

If no additional tape input, remaining card input as usual. 
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APPENDIX I ~ MULTIPOIMJLXONJLTRA INTS, 


EXAMPLE 1 



Node I and 2 are fixed. 

There is an applied displacement at node 3 in the negative y direc- 
tion. The distance between nodes 3 and 4 remains constant. 


Boundary Conditions 

011000000000 
000000000000 
021000000000 
Blank card 

Dependencies 

Node Dof 

3 2 

4 2 

4 3 


1-4 default condition 

12 4 
3 


Node • Dof 

3 2 

3 2 

3 3 


-.050 (specified displacement) 

1.0 restraint condition that 

1.0 3-4 remain rigid 


Note: In the boundary condition cards for node 4, degrees of freedom 2 
and 3 must be specified 0 because they are effectively eliminated from 
the solution (they are dependent degrees of freedom) . 
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EXAMPLE 2 


Rotation of Displacements at a Boundary 



For this problem we set the displacement u equal to zero and 
v free. This corresponds to the normal displacement along x 
equal to zero and the tangential displacement (in y direction) 
free. The global coordinates of the problem are- x, y. Thus: 

u = u cos q> - v sin f 

v - u s.in <p + v cos <p 

The multipoint constraint cards are: 

25 1 25 1 .86602540E00 25 2 -.50E00 

25 2 25 1 . 50E00 25 2 .8660254E00 

Note: On the boundary condition cards, since u = 0.0 node 25 

degree of freedom 1 is specified to be 0 but node 25 degree of 
freedom 2 is specified 1 or free since it remains an independent 
degree of freedom (although rotated) . If the normal displacement 
u were not fixed (and the transformation performed merely to apply 
a normal ' load) , it would have a 1 boundary condition specified. 
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APPENDIX II 


SOLUTION PACKAGE 


PODSYM - Solution of Symmetric Positive Definite Banded 
Matrix Equations 



Problem 

Solve AX = Y, where 

A = 


is a banded positive definite symmetric matrix; X is the desired 
solution vector; and Y is the known right side (load vector) . 
PODSYM is the user interface and supervisory routine. It uses the 
Cholesky algorithm which factors the total stiffness matrix into 
LL (where L is a lower priangular matrix) and then solves a 
pair of triangular sets of equations. 
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The factorization is carried out by subroutine QFACT which 
supervised the storage and data set allocation and subroutine 
QCHOL which generates the lower triangular L * matrix. QCHOL 
implements the Cholesky algorithm 'to factor the positive definite 
symmetric A matrix as the product of a lower triangular matrix 
with its transpose: 

T 

A = LL 1 

A straightforward argument establishes the possibility of decom- 
posing any positive definite matrix in this fashion. Once L has 
been obtained, it is not difficult to solve the system of linear 
equations AX = Y by calculating Z as the solution of the lower 

triangular system LZ = Y, and then X as the solution of the 

T 

upper triangular system L X - Z. The forward solution (LZ = Y) 

is accomplished by subroutines QFSOL and QFOR and the back 
T 

solution (L X = Z) by subroutines QBSOL and QBAC. Before the 

call to QBSOL, subroutine REVERS is called, which reverses the 

rows of L so that the last row becomes the first row, etc. This 

T 

is accomplished in order to sequentially access L during the 
back solution. 

The above algorithm is noteworthy for its assured stability 
and general efficiency. It is possible to carry out an error 
analysis of the procedure as it is represented on a digital compu- 
ter; such analysis shows that the computed L matrix satisfies 
an equation of the form 

A + E = LL T 

with bounds on the elements of E which show that E is small 

compared to A. The effect of rounding errors made in the subse- 

T 

quent solution of LZ = Y and L S = Z may then be taken into 
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account by (implicitly) introducing an additional perturbation F 
into A, and it is then concluded that the computed solution 
exactly satisfies the equation 


(A = E + F)X = Y 

Since E + F can be shown to be small, one would like to infer 
that X Q is almost an exact solution of the original equations, and 
unless A is too nearly singular, such a conclusion is indeed 
warranted. But if A is very ill-conditioned, no such result can 
be guaranteed, and X q may be far from the mathematically correct 
solution; in this event single-precision computation will not suf- 
fice for the calculation of an accurate solution, and since the 
solution will be very sensitive to small errors in A, it is un- 
likely that even a high-precision computation will yield satisfac- 
tory results unless A and Y are known (and supplied) to more 
than single -precis ion accuracy. The PODSYM subroutines make a 
fairly realistic attempt to detect and report pathological condi- 
tions of this sort. 

The large positive definite matrices that occur in practical 
work very often contain a large number of zero entries and the 
program seeks to benefit from the presence of these elements by 
modifying the standard Cholesky formulae. 


k-1 
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to read instead 


and 


k-1 


'kk 


\k- X 


'kj 


j— v (k) 


A., = 

xk 


A., 

xk 


k-1 

X 

j=ix(i,k) 


A.. I, . 
iJ k 3 


/ £ 


kk 
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APPENDIX III — CHANGING THE NUMBER OF NODES 
AND MEMBERS IN BEND 


The number of nodes and members allowed for a BEND problem 
is variable. The following changes must be made in the MAIN pro- 
gram for a successful up -dimensioning or down- dimensioning of nodes 
and/ or members . 

6 Reset MXMEM equal to the maximum number 
of members desired 

® Reset MXNOD to the maximum number of 
nodes desired 

9 Redimension C the main working area so 
that it is greater than the greater of 
10.5 * MXNOD + 6.5 * MXMEM or 24 * MXNOD 
+ 1975 

® Set NCORE equal to the dimension of the 
C array 

• For maximum efficiency set MHCON equal 
to a prime number approximately 1/3 more 
than MXMEM 

9 For maximum efficiency set NHCON equal to 
a prime number approximately 1/3 more than 
MXNOD. A table of prime numbers appears in 
Ref. 6, page 870. 

® For CDC versions, dimension IEAC (random 
access array) toMXMEM+1. 

The following limitations on the number of property specifica- 
tions exist in BEND, and can only be changed through redimensioning 
within the program: 
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Current Limitation 


Property or Section 

Number of MATl and MAT2 property 20 

types 

Number of MBM property types 20 

Number of beam section property 100 

types 

Number of types of edge loads 100 

for EDG1, EDG2j EDG3 arrays 

Number of types of surface 100 

loads applied (SURF) 

Number of multipoint 500 

constraint conditions 

Number of beam element 100 

distributed load specification types 

Number of boundary condition 100 

type specifications 
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SYMMETRY 


APPENDIX IV EXAMPLES OF PLATE BOUNDARY CONDITIONS 



BOUNDARY CONDITION: 1 - FREE; 0 = FIXED 


NODE 

U 

V 

w 


H 


8 * 

e v 

7 

xy 

*XX 

K vy 

K xy 

1 

0 

0 

1 



0 

1 

1 

0 

1 

1 

0 

2 

1 

0 

1 

0 

1 

0 

1 

1 

0 

1 

1 

0 

3 

1 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

0 

4 

0 


1 

1 

0 

0 

1 

1 

0 

■ 1 

1 

0 

6 

1 

1 

0 

0 

1 

o' 11 

1 

1 

o (2) 

0 

0 : 

1 

6 

1 

1 

0 

0 

1 

1 

1 

1 

1 

0 

0 • 

1 

7 

0 

1 

0 

1 

0 

0 

1 

1 

0 

0 

0 

0 

8 

1 

. 1 

0 

i 

1 

0 

o (,i 

1 

1 

0 I2) 

0 

0- 

1 

8 

1 

1 

0 

1 

0 

1 

1 

1 

1 

0 

0 

1 

9 

1 

1 

0 

0 

0 

0 U) 

1 

1 

0<2) 

0 

0 

1 

9 

1 

1 

0 

0 

0 

* 7 

1 

1 

1 

0 

0 

1 


ALL ZERO ALL ZERO IF 

•IF NO IN-PLANE LOAD NO IN-PLANE LOAD 

NOTES- 

{II = 0 FOR ZERO OR UNIFORM DIRECT IN-PLANE LOADS 
« 1 FOR NON-UNIFORM IN-PLANE LOADS 
(21 = 0 FOR ZERO IN-PLANE SHEAR LOAD 

= 1 FOR. IN-PLANE SHEAR 
NODE 5 COMPLETELY FREE 
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BOUNDARY CONDITION: 1 


NODE 

U 

V 

w 

0 

X 

0 

V . 

0 

z 


1 

0 

0 

1 

0 

0 

0 


2 

1 

0 

1 

o' 

1 

0 


3 

0 

0 

0 

0 

0 

0 


4 

0 

1 

1 

1 

0 

0 


6 

0 

0 

0 

0 

0 

0 


7 

0 

0 

0 

0 

0 

0 


8 

0 

0 

0 

0 

0 

0 

: 

9 

0 

0 

0 

0 

0 

0 



NODE 5 COMPLETELY FREE 
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o CLAMPED - EDGES RESTRAINED 
AGAINST MOVEMENT IN THE X-Y 
PLANE 


LATERAL 



SECTION A-A 


A 


SECTION B-B 


£ ^ X, U SYMMETRY 


- FREE; 0 = FIXED 


'xy 


VV 


xv 


_0 

0 


































ALL ZERO IF ALL ZERO IF 

MO INPLANE LOAD NO INPLANE LOAD 


NOTES 

(II = 0 FOR ZERO OR UNIFORM DIRECT IN-PLANE LOADS 
= 1 FOR NON-UNIFORM IN-PLANE LOADS ' 

(2) = 0 FOR ZERO IN-PLANE SHEAR LOAD 
= 1 FOR IN-PLANE SHEAR 
NODE 5 COMPLETELY FREE 
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• SIMPLY - SUPPORTED: EDGES RESTRAINED 






SECTION A-A 


SECTION 3-B 


1 


x,u SYMMETRY 


BOUNDARY CONDITION: 1 - FREE; 0 = FIXED 


NODE 

U 

V 

w 

°x 

°v 

0 

z 

e x 

e v 

^xv 

*xx 

\v 

*xv 

1 

0 

0 

1 

0 

0 

0 

1 

1 

0 

1 

1 

0 

2 

0 

0 

1 

0 

1 

0 

1 

1 

0 

1 

1 

0 

3 

0 

0 

0 

0 

1 

0 

1 

,11) 

0 , 

0 

0 

0 

3 

0 

0 

0 

0 

1 

0 

1 ; 

0 (2) 

0 

0 

0 

0 

4 

0 

1 

1 

1 

0 

0 

1 

1 

1 

1 

0 

1 

1 

0 

6 

1 

id) 

0 

0 

1 ! 

0 

1 

1 (1 > 

o (1) 

0 ’ 

0 

1 

6 

1 

0 (2) 

0 

0 

1 

' 0 

1 

0 (2) 

1 ( 21 

0 

0 

1 

7 

0 

0 

0 

1 

0 

0 

id) 

1 

0 

0 

0 

0 

7 

0 

0 

0 

1 

0 

0 

0 ,2 » 

L 1 

0 

0 

0 

0 

8 

1 (1 > 

0 

0 

1 

o 

0 

i 1 {1 > 

1 

0 (1) 

0 

0 

1 

8 

ESI 

0 


1 

0 

0 

0<2) 

1 

1 (2! 

0 

0 

1 

9 

0 

0 

0 

0 

0 

Q 

,d> 

1*1» 

0 

0 

0 

1 

9 

0 

0 

0 

0 

0 

0 

0 (2> 

o <2 > 

0 

0 

0 

1 


NOTES- 

{1} WHEN DISPLACEMENTS ARE NOT RESTRAINED IN THE DIRECTION OF THE EDGE 
(21 WHEN DISPLACEMENTS ARE FULLY RESTRAINED 
NODE 6 COMPLETELY FREE 
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SECTION 3 


INSTRUCTIONS FOR USE OF HEX 

A Program for the Elastic-Plastic Analysis 
of Three Dimensional Solids 



INSTRUCTIONS FOR THE USE OF HEX 


HEX is a finite element program to treat the elastic, elastic- 
plastic or elastic-cyclic plastic response of arbitrary three dimen- 
sional solid structures. The program uses a family of isoparametric 
hexahedra elements (Refs. 1 and 2) consisting of eight-node hexahedra 
and hexahedra with up to 12 additional midside nodes as shown in 
Fig. 1. Sample problems from the present program can be found in Ref. 3. 

The program is capable of treating the elastic and the elastic- 
ideally plastic response of orthotropic materials. In addition, 
consideration is given to isotropic materials exhibiting elastic- 
ideally plastic, linear strain hardening, or nonlinear strain 
hardening behavior- Further, the kinematic hardening theory of 
plasticity is used (Refs. 4-6) so that provision for cyclic loading 
conditions involving reversed plastic deformation is included. 

The input to the program is categorized in the following 
sections : 

I. Problem Title FORMAT (20A4) 

Any 80-character title describing the problem. 

II. NPNTC, NPRNT, IRESRT , NUTAP, INPRT FORMAT (515) 

0 < NPNTC <63: 

NPNTC is the sum of the following integers corresponding to 
the option desired. 

If NPNTC = 0 No intermediate printout 

= L Print the load vector 
= 2 Print element stiffness matrix 
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X 

Local Coordinates Cartesian Coordinates 

Fig. 1 Isoparametric Hexahedra 

NPNTC = 16 Print each element stiffness matrix entry 
to be stacked with its stacking index 
= 32 Print the total stiffness matrix 

For example, if it is desired to print the load vector and 
the total stiffness matrix, NPNTC - 1+ 32 = 33. 

NPRNT: 

If <0, perform elastic analysis pnly. 

If >0, perform plastic analysis, printing output every 
NPRNT increments of load. 
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IRESRT: (See Section XVIII) 

0, Elastic (and/or) plastic run. Do not generate 
a restart tape 

1, Elastic (and/or) plastic (cyclic) run. Build a 
new restart tape 

2, Plastic run with elastic values from previously 
created restart tape 

3, Plastic restart run starting at some specified 
load level 

4, Cyclic restart run 
NUTAP: 


Applicable only if restarting from a restart tape • 
0, No new tape written 


1, Complete new tape created and additional restart 
data written 


INPRT : 

O', Write restart tape only at P = PMAX (i.e., at maximum 
load) 

N, Write restart tape every N increments of load 


III. Node Specification (1^15) 


This section defines an allowable set of external node point 
numbers. The maximum node number that can be used is 9999, The 
program uses this information in two ways. First to set up a table 
of allowable node points that is used to check all subsequent node 
point input. Secondly, the program converts each external node num- 
ber to an internal number consecutively in the order that the node 
appears on the input card. Consequently the order of the input of 
external node numbers is completely arbitrary and need not be in- 
creasing monotonically . In practice the node numbers should be num- 
bered so as to minimize the bandwidth. Once the input is read the 
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program operates with the internal numbers which are now numbered 
from 1 through the number of nodes in the model. In this manner 
the node ordering and therefore the bandwidth of the stiffness 
matrix can be easily changed and nodes can be inserted or deleted 
by changing the external node specification. 

The input is specified by entering the appropriate number on 
the input cards in fields of five. However, for this purpose the 
user can also utilize a shorthand form of the input. That is, spe- 
cifying m and -n consecutively is the equivalent of the speci- 
fication of nodes m, m + 1, m + 2, ..., n and specifying m, -p, 
and -n consecutively is the equivalent of the specific tion of 
nodes m, m + p, m + 2p, ..., m + kp where m + kp is the highest 
integer of the form less than or equal to n. For example,, the spe- 
cification of nodes 1 through 100 is written as 1 - 100 and 
1, 3, 5, . .., 99 as 1-2-99. This card input appears in fields of 
5 (15 Format) with 16 items per card. Any number of continuation 
cards may be used. A zero or blank 15 field ends the specification . 

rv ’ Member Connectivity (Node Nu mbers of Each Member) 

FORMAT (A4, 6X, 915) ~ ~ ~ 

The first alphanumeric field defines the element type as 1 fol- 
lows 

HEX8 - eight node hexafyedra element 

HX20 - hexahedra element with up to twelve additional 
midside nodes 

MSGN - mesh generator for eight-node hexahedra 

SEND - ends the section 

The first integer field designates the member number. The next 
integer field designates the connecting nodes as follows; 

HEX8 - The eight corner nodes are specified in the 
order shown on Fig . 2 . 
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HX20 


Note: 


MSGN 


- The eight corner nodes are specified first 

in the order shown in Fig. 2. The additional 
midside nodes are specified on a second card 
(FORMAT (1215)) according to the order shown 
in Fig.. 3. An element can be constructed 
with any combination of midside nodes. If a 
midside node is not specified, a zero is 

placed in its corresponding position on the 
card. If a midside node is specified for 
an element, a midside node should also be 
specified for every element coincident with 
that side . 

The hexahedra elements with midside nodes 
can have curved boundaries . The shape of 
any side is described by a quadratic poly- 
nomial through the three node points . If 
a side containing a midside node is 
straight, the coordinates of the midside 
node - can be read in as zero. In this 
case the nodal coordinates are calculated 
within the program. 

- The nodal connectivity is automatically 
generated for eight-node hexahedra accord- 
ing to the. convention in Fig. 2. The 
first card for this input option requires 
the four -character identifier MSGN. A 
second card is. required with FORMAT (515) 
which specifies the number of element 
subdivisions in three coordinate direc- 
tions. The integers in the five fields are 
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the subdivisions in the x^, x^, direc- 

tions, and the first node number and element 
number in the mesh. The coordinate system 
must be right handed. The mesh generator 
assumes that the nodes are numbered consecu- 
tively along the x^ direction specified 
in the first 15 field then incremented on 
the x^ - x^ face until that face is com- 
pletely specified. • The mesh is then incre- 
mented in the direction and then con- 

secutively numbered on the x^ - x^ plane. 

The element numbers are numbered consecutively 
starting at the specified starting value in 
the positive x^ direction. The mesh gener- 
ating capability still requires that the nodes 
be defined in Section IV and that the nodal co- 
ordinates be explicitly defined in Section VI . 

An example of the convention for a 4x3x2 
mesh in Cartesian coordinates with the starting 
node and element equal to one is shown in Fig. 4. 

V . X -Coordinates of Nodes FORMAT (E15.7, 1315) 

The global x-coordinates of the nodes appearing in the 15 
fields are set to the value in the E15.7 field. Any number of 
continuation cards may be used; their first fifteen columns are 
ignored. A zero or blank 15 field terminates the card scan for 
a given x-coordinate . A zero or blank first 15 field (columns 
16-20) on a noncontinuation card ends the section. Both shorthand 
representations of Section III are allowed. 

VI. Y-Coordinates of the Nodes . Same as Section V. 

VII. Z-Coordinates of the Nodes. Same as Section V. 
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Fig. 4 Node and Member Ordering Convention 
for Automatic Mesh Generation 


VIII. Boundary Conditions FORMAT (13, 211, 5X, 1415) 

The first three fields give the boundary conditions specifica- 
tions in the order: u, v, w; where u, v ; w are global displace- 
ments in the global x, y, z directions, respectively. 

Zero denotes a fixed degree of freedom 

One denotes a free degree of freedom 

Two will result in the application of a unit general- 
ized displacement, or a corresponding card may be 
included in Section IX (dependent degrees of free- 
dom) designating the magnitude of the applied gen- 
eralized displacement. 

The 1415 fields give the applicable nodes for the boundary 
condition specifications, with both shorthand notations of Sec- 
tion III permitted. Any number of continuation cards may be used 
for a given specification. However, only the 1415 fields are 
used on a continuation card. A zero or blank 15 field terminates 
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the card scan for a given boundary condition specification . Note: if 
the last field of a card (columns 7^-80) is the last specification, 
an additional blank card (continuation card) must follow. A zero or 
blank first 15 field (columns 11-15) on a noncontinuation card ends 
the section. If a node’s boundary conditions are not specified in 
this section, all the degrees of freedom are assumed to be free. To 
change this default condition, the first card of this section should 
be set to the desired default with all nodes used in the problem 
specified. Note: Maximum number of nodes is currently 2500. 


IX. Dependent Degrees of Freedom FORMAT (215, 2(5X, 215, E15.7)) 

This section designates the input for both single and multipoint 
constraints as well as applied displacement of the form: 


1 ) 


F 'i • V'l + a 2 5 2 + 


a b 
n n 


where b^ is a dependent degree of freedom, 

f>l • • • f are independent degrees of freedom, 

and a,, a„, •**, a are coefficients. 

I 5 2 n 


2) Rota tion of displacements at a node 


b . = a, b . + a„ b . + a 0 5, 
x 1 r Z j ok 


*\i dy 'P 

8 j " P l 6 i + P 2 6 j + P 3 S k 


h = y ft. + 7 b. + ‘y„h 

k 1 l 2 j 3 k 


where the a, (3 , 7 ‘s are the direction cosines of 

the rotation, 6., b, are the displacements with 

respect to the original global directions and 5., 
r \/ 

fh , are the components of displacements at the 

node with respect to the new coordinate axes . An 
example of this capability is given in Appendix I. 
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3) Applied generalized displacement 



where the coefficient is the applied generalized 

displacement. 

The first two 15 fields designate a node number and a degree of 
freedom (i.e., 1-3). The dependency is defined in the following 
three fields. The two integer fields designate the node number 
and degree of freedom number and the coefficient by the floating 
point field. If there is another dependency for the node, it is 
designated in a similar fashion in the next three fields. Any 
number of continuation cards can be used with the first two fields 
blank. The section is ended by a blank or zero in the third in- 
teger field (blank card) . Examples of the use of multipoint con- 
straints are in Appendix I. 

X. Material Properties 

The input for material properties consists of cards as fol- 
lows : 

CARD 1: FORMAT (A4, IX, 5E15.7) 

MATl ~ four-character identifier 

EONE = Young's modulus in principal property axis (1) 

ETWO = Young ' s modulus in principal property axis (2) 

ETEE = Young's modulus in principal property axis (3) 

VONTO = Poisson's ratio, 

VTOTE = Poisson's ratio, v ^ 

CARD 2: FORMAT (4E15.7) 


VTETO 

= Poisson's ratio. 

V 32 

GONTO 

= Shear 

modulus , 

G 12 

GTOTE 

= Shear 
1 

modulus , 

G 23 

GONTE 

= Shear 

modulus , 

G 13 
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CARD 3: FORMAT (3E15.7) 


BETA1 = I Angles between 1, 2, 3 material axes 

BETA2 = V and local element x, y, z axes (defined 

BETA3 = ] in the note at the end of the section) . 

CARD 4: FORMAT (3E15.7) 

* 

TALF-1 Coefficient of thermal expansion in 

the principal 1 direction 
TALF-2 Coefficient of thermal expansion in 

the principal 2 direction 
TALF-3 ' Coefficient of thermal expansion in 

the principal 3 direction 

CARD 5: FORMAT (5E15.7) 

SIGOX = Yield stress in principal 1 direction 

SIGOY — Yield stress in principal 2 direction 

SIGOZ = Yield stress in principal 3 direction 

SIGOXY = Shear yield stress in principal 

1- 2 plane 

SIGOYZ = Shear yield stress in principal 

2- 3 plane 

CARD 6: FORMAT (5E15.7) 

SIGOZX = Shear yield stress in principal 

3- 1 direction 

RMOSN = .If RMOSS + 0; RMOSN = n, the shape 

parameter used in Ramberg -Osgood repre- 
sentation of stress -strain behavior. 

If RMOSS = 0; RMOSN = a the slope of 
the linear strain -hardening stress-strain 
representation, i.e., a = E^,/E where E^ 
is the tangent modulus . 
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RMOSS 


If RMOSS ^ 0; RMOSS = Ramberg -Osgood 
parameter ^ 

Note (1) : If RMOSN - 0 and RMOSS = 0 

the material for the element (s) is assumed 
to be elastic-ideally plastic. 

YLDST = Yield stress (for isotropic materials) 

Note (2): Only initially isotropic mate- 

rials can be treated when considering 
linear or nonlinear strain- hardening. 


RMOSE = Ramberg-Osgood parameter E (Young’s modulus) 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 


each member list ins 


Note: The transformation between the 

1, 2, 3 material axes and the local ele- 
ment x, y, z axes are defined by path 
dependent angular rotations P 3 • 

First a rotation about the 1 axis, 
then a rotation about the 2 axis in 
the transformed system and then a rotation, 
f'o about the 3 axis in the transformed 
system 
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ACUM 


Table describing yield stress versus accumulated plastic 
strain relationship 

FORMAT (A4/A4, IX, 15, /, A4, IX, 2E15.7, /, (1615)) 

CARD 1: 


ACUM - four character identifier 

CARD 2: 

ACUM - four character identifier 

NUMP ~ number of pairs of points used to repre- 
sent the yield stress versus accumulated 
plastic strain 

CARD 3: 


ACUM - four character identifier 

VARS'iG - yield stress 

SUMEPS - accumulated plastic strain 


Card 3 is repeated "NUMP” times where NUMP is the number of pairs 
in the table. A maximum of 20 pairs is allowed. 

Succeeding cards give applicable members; both shorthands of 
Section III are permitted. Any number of continuation cards may 
be used for a given specification. A zero or blank 15 field ends 
each member listing. 
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GAUSS or LOBATTO - Designation of stress points within an element 


CARD 1: FORMAT (A4, IX, 3H5) 

C GAUSS = five letter character identifier 


or 


LOBAT 

NX 


five letter character identifier 

number of stress points in local x 
direction 


NY = number of stress points in local y 

direction 

NZ = number of stress points in local z 
direction 


CARD 2: Applicable members; both shorthands of 

Section III are permitted. Any number of con- 
tinuation cards may be used for a given spe- 
cification. A zero or blank 15 field ends 
each member listing. 

SEND in columns 1-4 ends the section. 


Note: Up to 8 Gauss or 10 Lobatto points are per- 

mitted in any one direction and a total of 
50 points are permitted in any one element. 
For local coordinate definitions see Figure 2. 
For a definition of these integration points 
see Ref. 7. 
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XI. Applied Loads 

Three different types of loading are currently admissible — 
concentrated loads at nodes, distributed loads on a surface of a 
hexahedron member, and thermal loads. Each loading situation is 
designated by one of the following four-character key words: 

CONC (Concentrated load) 

SURF (Surface distributed load) 

TEMP (Temperatures for thermal loading) 


Distributed load - 


CARD 1: FORMAT (A4) 

Key word SURF 

CARD 2: FORMAT (15, 5X, 815) 

First 15 field: number of nodes on the loaded 

surface. There is a minimum of 
four nodes per surface and a 
maximum of eight in order to 
accommodate midside nodes. 


Second through' 

fifth 15 fields: These four fields contain the 

four corner node numbers in 
counterclockwise or clockwise 
order. 
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Sixth through 

ninth 15 fields: These four fields contain the 

midside node numbers (if any) . 
If a particular midside node 
is omitted, a zero or blank is 
placed in the appropriate field 
on the input card . 

The input order for a given surface is shown below for counter- 



The magnitude of_the surface tractions at each node specified above 
follow: There is one card for each node specified. 

SUCCEEDING CARDS: FORMAT (15, 3E15.7) 

NUM = Node number specified on CARD 2 (15 field) 

TX = Surface traction in global x-direction 

TY = Surface traction in global ydirection 

TZ = Surface traction in global z-direction 

A blank card ends the input for surface loads. 

Note; The keyword SURF appears only once at the beginning 

of the specification of surface loads. Additional 
input for surface loads begin with card 2. 
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Concentrated load 


CARD 1: FORMAT (A4) 

Key word CONC 

CARD 2: FORMAT (15, 3E15.7) 

NUM = Node number at point of application of 
load (15 field) 

FX = Load component in the global x-direction 
(first E15.7- field) 

FY = Load component in the global y-direction 
(second E15 . 7 field) 

FZ = Load component in the global z -direction 
(third E15.7 field) 

As many load components as desired may be specified. A' zero or 
blank first 15 field ends the specification of concentrated 
loads . 

Note: The keyword CONC appears only once at the beginning of the 
specification of concentrated loads. Additional input for 
concentrated loads begin with card 2. 

Thermal load - 

Each card after the first designates the temperature increment 
from some reference temperature followed by its applicable nodes. 
The input format is- identical to that described for the input of 
nodal coordinates in Section V. 

CARD 1; FORMAT (A4, IX, E15.7) 

Key word TEMP 

TREF = Reference temperature in E15.7 field 

The succeeding cards are as in Section V for nodal coordinates. 
That is, the nodal temperature in the first E15.7 field followed 
by the applicable nodes. A blank card ends the specification of 
nodal temperatures . 
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Note: The nodal temperature specified at the corner nodes of an 

individual hexahedra are interpolated to obtain the temperature 
used in calculating stresses and strains at the stress points. 

SEND in columns 1-4 ends the section for applied loading. 

XII . Members to be Printed - FORMAT (1615) 

Specify the members whose strains and stresses are to be printed. 
Both shorthands of Section III are allowed. A maximum of 2500 mem- 
bers may be specified. Members in excess of 2500 and undefined mem- 
ber numbers are ignored. A blank card or card with only zero entries 
ends the section. 

XIII . Nodes to be Printed - FORMAT (1615) 

Up to 2500 nodes whose displacements are to be printed for the 
analysis, as per Section XII. A blank card ends the section. 

XIV. Parameters for Plastic Analysis - FORMAT (2E15.7) 

PMAX - Maximum load to be applied for this half cycle 

PPCT = Load increment as percentage of yield load 

XV. Parameters for Succeeding Load Cycles - FORMAT (A5, 15, 2E15.7) 

CARD 1: FORMAT (A5) 

Keyword defining the load type as follows: 

CYCLE - Same load distribution as last half cycle 
NEW *= New load distribution to be specified 
If any other word (or a blank card) is used the analysis is complete. 

CARD 2: FORMAT (15, 2E15.7) 

NPRNT = If equal to zero, no additional load cycle and end 
of problem. If greater than zero, print output 
every NPRNT increments. 

PMAX = New maximum load for this load cycle 
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Note : 


PPCT = Load increment as percentage of yield load. 

When CYCLE is specified, the original load distribution is 
maintained. In this case PPCT should be set equal to one- 
half the original value since the elastic range for subse- 
quent cycles is twice the elastic range for initial loading. 

XVa. If NEW is used in Section XV a new load must be specified. 
This is accomplished with input Identical to Section XI - 
Applied Loads. If CYCLE is specified in Section XV this 
section is omitted. 


XVI. Change Plastic Material Properties When Cycle is Specified 
CARD 1: FORMAT (A4, IX, 15) 

MATl = Key word 

IMEM = One member number from the group that is 
to be changed. 

CARD 2: FORMAT (5E15.7 


SIGOX 
SIGOY 
SIGOZ \ 

SIGOXY 
SIG0YZ - 

CARD 3: FORMAT (5E15.7) 

SIGOZX 
RMOSN 
RMOSS V 
YLDST 
RMOSE j 


Defined in Section X 


Defined in Section X 


The section is ended with a section end card, SEND. 
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Note: If there is no change in material properties, a SEND- card 

must appear in the appropriate place in the input stream. Suc- 
ceeding load reversals (CYCLE) or new loads (NEW) are specified 
with additional input beginning with Section XV. A blank card 
(or any alphnumeric other than CYCLE or NEW) branches to Sec- 
tion XVII to end the job. 


XVII. 


Each problem's input must be ended with a card reading "ENDb" 
where "b" denotes a blank, in columns 1-4. The last problem in 
a run should end with a card reading "STOP" in column 1-4 instead 
of "ENDb. 11 


XVIII. Restarting a Problem 

The initial restart tape is created on Unit 21. Subsequent 
restart jobs mount the restart tape as Unit 21 and if desired 
(NUTAP > 0) copies and continues the restart tape on Unit 22. 


CARD 1: Problem Title - FORMAT (20A4) 

As i n Se ction I 

CARD 2: NPNTC , NPRNT, IRESRT , NUTAP, INPRT - FORMAT (515) 

As in Section II. Here NPNTC is ignored. 

If IRESRT is equal to 2 

Restart from an elastic critical load using a 
previously created restart tape 

CARD 3: PMAX, PPCT - FORMAT (2E15. 7) 

PMAX = Maximum load to be applied for this half cycle 

PPCT = Load increment as percentage of yield load 

Input is then continued from Section XV 
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If IRESRT is equal to 3 

Continuation from some previously plastic load. The 
job can be restarted from a previous maximum load or 
from any intermediate load level. 

CARD 3: PMAX, NRSRT - FORMAT (E15.0, 15) 

If starting from a previous maximum load NRSRT can 
be ignored. If starting from some intermediate 
load level, NRSRT is a unique number obtained as 
output from the job that generated the restart tape. 
This output is of the form, 

k -k •k k >V 

PLASTIC ANALYSIS VALUES FOR RESTART TAPE RAVE BEEN 
WRITTEN AT P - 1. 5000000E+03 FOR PRINT INTERVAL 30 

* -jV * * * 

In this case NRSRT is equal to 30 and the load is 
not necessarily the last maximum value. 

Note: For cases when the job ends abnormally (for example when 

time or lines are exceeded) the job can be restarted v?ith 
a value for NRSRT obtained from the output. In this case 
the last value of NRSRT should not be used. This will 
insure that a complete set of data will be accessed from 
the restart tape. 

If IRESRT is equal to 4 

Restart from a previous maximum load level followed 
by application of a new load (NEW) or reversal of a 
previous load distribution (CYCLE). A new load or' 
a reversal of a previous load distribution can be 
placed on the structure from any previously generated 
and saved state of stress. 
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CARD 3: 


PMAX, NRSRT - FORMAT (E15.0, 15) 

If starting from a previous maximum load NRSRT can 
be ignored. If starting from some intermediate load 
level, NRSTR is a unique number obtained as output 
from the job that generated the restart tape. This 
output has the form shown in the preceding dis- 
cussion for IRESRT equal to three. The value of 
PMAX is ignored. 

Succeeding cards follow as in Sections XV through XVII. 

XIX. Alternate Input Method 

All or part of the input may be read from a tape. This tape is 
mounted as Unit 23 with the appropriate job control cards. The tape 
is accessed by specifying the keyword, TAPE (FORMAT A4), in the proper 
place in the- card input stream. Input will then be read from the tape 
until a SEND (A4) is encountered on the tape. This returns the input 
reader to the card input stream. 

Example.: 

Input for member connectivity. Section IV, and coordinates have 
been generated on a tape in the prescribed manner. 

Input cards are as follows: 

TITLE 

0 5 

1-100 

TAPE (reads data from a tape) 
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HEX8 

. (additional elements if desired) 

* 

SEND 

TAPE (reads X-coordinates from tape) 

(Blank card to end section for X-coordinates) 
TAPE (reads Y-coordinates from tape) 

(Blank card to end section for Y-coordinates) 
5.0 (card input for Z -coordinates 

* 

(Blank card ends section for Z-coordinates) 
If no additional tape input, remaining card input as usual. 

Note: All SEND formats used to end sections must be in 

card form as the SEND encountered on the tape only sends the 
reader back to the card input stream. Also blank cards which end 
sections are best not put on the tape but rather on cards as this 
allows the user to add card input to tape input, i.e. , if a blank 
ends a section and is encountered on the tape format then the sec- 
tion ends and no new data may be added, but if the blank does not 
occur on tape the section is still open and can be closed with a 
blank card in the card input stream. 
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APPENDIX I — MULTIPOINT CONSTRAINTS 



Node 1 and 2 are fixed. 

There is an applied displacement at node 3 in the negative y direc- 
tion. The distance between nodes 3 and 4 remains constant. 

Boundary Conditions 

011000 1-4 default condition 

001000 124 

021000 3 

Blank card 

Dependencies 


Node 

Dof 

Node 

Dof 



3 

2 

3 

2 

-.050 

(specified displacement) 

4 

2 

3 

2 

1.0 

restraint condition that 

4 

3 

3 

3 

1.0 

3-4 remain rigid 


Note: In the boundary condition cards for node k, degrees of freedom ? 

and 3 must be specified 0 because they are effectively eliminated 
from the solution (they are dependent degrees of freedom) . 



EXAMPLE 2 Rotation of Displacements at a Boundary 



For this problem we set the displacement u equal to zero and 
v free . This corresponds to the normal displacement along x 
equal to zero and the tangential displacement (in y direction) 
free. The global coordinates of the problem are x, y. Thus: 

u = u cos <p - v sin qp 

v = u sin 9 + v cos <p 

The multipoint constraint cards are: 

25 1 25 1 .86602540E00 25 2 -.50E00 

25 2 25 1 . 5OE00 25 2 .8660254E00 

Note: On the boundary condition cards * since u = 0.0 node 25 

degree 'of freedom 1 is specified to be 0 but node 25 degree of 
freedom 2 is specified 1 or free since it remains an independent 
degree of freedom (although rotated) . If the normal displacement 
u were not fixed (and the transformation performed merely to apply 
a normal load ) , it would have a 1 boundary condition specified. 
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APPENDIX II 


SOLUTION PACKAGE 


PODSYM - Solution of Symmetric Positive Definite Banded 
Matrix Equations 


| PODSYM] 


I 1 

i QFACT) | Qf'solW ^OFOR 'I 

f • 

l 

rO GHOL I ■ REVERS j 

\ 

IqbsolW — «» fQBAC 1 


{■RE VERS i 

Problem 

Solve AX = Y, where 


A = 



is a banded positive definite symmetric matrix; X is the desired 
solution vector; and Y is the known right side (load vector) , 
PODSYM is the user interface and supervisory routine. It uses the 
Cholesky algorithm which factors the total stiffness matrix into 
LL (where L is a lower triangular matrix) and then solves a 
pair of triangular sets of equations. 
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The factorization is carried out by subroutine QFACT which 
supervised the storage and data set allocation and subroutine 
QCHOL which generates the lower triangular L matrix. QCHOL im- 
plements the Cholesky algorithm to factor the positive definite sym- 
metric A matrix as the product of a lower triangular matrix with 
its transpose: 

T 

A = LL 

A straightforward argument establishes the possibility of decompos- 
ing any positive definite matrix in this fashion. Once L has been 
obtained, it is not difficult to solve the system of linear equa- 
tions AX = Y by calculating Z as the solution of the lower tri- 
angular system LZ = Y, and then X as the solution of the upper 

T 

triangular system L X = Z. The forward solution (LZ = Y) is 
accomplished by subroutines QFSOL and QFOR and the back solu- 
tion (L^X = Z) by subroutines QBSOL and QBAC. Before the call 
to QBSOL, subroutine REVERS is called, which reverses the rows 

of L so that the last row becomes the first row, etc. This is 

T 

accomplished in order to sequentially access L during the back 
solution. 

The above algorithm is noteworthy for its assured stability 
and general efficiency. It is possible to carry out an error 
analysis of the procedure as it is represented on a digital compu- 
ter; such analysis shows that the computed L matrix satisfies an 
equation of the form 

A + E = LL T 

with bounds on the elements of E which show that E is small 

compared to A. The effect of rounding errors made in the subse- 

T 

quent solution of LZ = Y and L X = Z may then be taken into 
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account by (implicitly) introducing an additional perturbation F 
into A, and it is then concluded that the computed solution X 

o 

exactly satisfies the equation 


(A = E + F)X = Y 

Since E + F can be shown to be small, one would like to infer 
that X q is almost an exact solution of the original equations, 
and unless A is too nearly singular, such a conclusion is indeed 
warranted. But if A is very ill-conditioned, no such result can 
be guaranteed, and X^ may be far from the mathematically correct 
solution; in this event single -precis ion computation will not suf- 
fice for the calculation of an accurate solution, and since the 
solution will be very sensitive to small errors in A, it is un- 
likely that even a high -precision computation will yield satisfac- 
tory results unless A and Y are knox\m (and supplied) .to more 
than single-precision accuracy. The PODSYM subroutines make a 
fairly realistic attempt to detect and report pathological condi- 
tions of this sort. 


The large positive definite matrices that occur in practical 
work very often contain a large number of zero entries and the 
program seeks to benefit from the presence of these elements by 
modifying the standard Cholesky formulae 


k-1 


'kk 


- 1 


’l 

j=l J 


and 


f 'ik 


k-1 

A., - ) l. . L . 

ik L 13 kj 

3 = 1 


for i > k 
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to read instead 


and 



\k 


k-l 


1 

j-V(k) 


i 

? 



k-l 

A ik “ 1 ^ij^kj 

j=M-(i,k) 
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SECTION 4 


INSTRUCTIONS FOR 
THE USE OF REVBY 


A Program for the Elastic-Plastic 
Analysis of Bodies of Revolution 



INSTRUCTIONS FOR THE USE OF REVBY 


REVBY is a finite element program to treat the elastic, 
elastic-plastic, or elastic -cyclic plastic response of ortho- 
tropic axisymmetric solids of revolution under axisymmetric 
loadings. Currently, three different types of elements are 
available for the analysis of thick and thin bodies of revolution: 


1 . Revolved Triangular Element 

The revolved triangular element used is similar to the 
element presented in Ref. 1, restricted to axisymmetric 
loading, i.e., independent of the circumferential direc- 
tion, but including orthotropic material properties (see 
Fig. 1). 

2 . Thin Shell Element 

The thin shell element used was developed at Grumman and 
the theoretical derivations are presented in Refs . 2 and 
3. Cubic Hermitian polynomials represent the three 
displacement components. Again, only axisymmetric load- 
ing is allowed but orthotropic material properties are 
included (see Fig. 2). 

3. Thin Ring Stiffener Element 

The thin ring stiffener element is similar to the one 
presented in Ref. 4 with the exception that the ring may 
be attached with arbitrary eccentricities . No warping 
of the cross section is allowed and the shear center and 
centroid must coincide. Five different cross-sectional 
geometries are available. 
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Displacement Assumption: 


u = 

a. 

■+■ a n r 

+ 

a„z 

r 

1 

2 

3 

u ■ = 

a. 

•+ -a c r 

+ 

a^z 

z 

4 

5 


6 

u _ - 

Cl 

a 7 

+ a g r 

4- 

a g z 


Initial Strain Distribution: 
e..(r,z) = constant 


Fig. 1 Axisymmetric Revolved Triangle 



Displacement Assumption: 


■ u = 

a i 

+ 

a 2^ 

+ 

. 2 
a 3^ 

4“ 

V 

VJ =s 

b i 

+ 

V 

+ 

, .2 
V' 

+ 

V 

V = 

c i 

+ 

c 2~~ 

+ 

-2 
C ^ *5 

+ 

> 

c c. 

c 4 


Initial Strain Distribution: 

'«i«> ' f4< £ > I 

{ 4co l + *5(0 < 

l -ie«> I 



Fig.- 2 Axisymmetric Thin Shell Element 



The plasticity theory used is presented in Refs. 5, 6, 
and 7. The input to the program is categorized in the fol- 
lowing sections . 

I. Problem Title (20A4) 

Any 80-character title describing the problem. 


II. NPNTC and NPRNT (215) 


NPNTC: is for selective printing. NPNTC is the sum of the 

following integers corresponding to the options desired. 


If NPNTC 


0, no intermediate output 

1, print total load vector 

8, print nonzero element stiffness matrix 
entries with their stacking indices 

16, print the total stiffness matrix 


NPRNT: If <. 0, perform elastic analysis only 

If > 0, perform elastic and plastic analysis, 
printing the results of the plastic 
analysis every NPRNT load increments 


HI. Nodes by Partition (1615) 

This section defines an allowable set of external node point 
numbers. The maximum node number that can be used is 9999. The 
program uses this information in two ways. First to set up a fable 
of allowable node points that is used to check all subsequent node 
point input. Secondly, the program converts each external node num- 
ber to an internal number consecutively in the order that the node 
appears on the input card. Consequently the order of the input of 
external node numbers is completely arbitrary and need not: be in- 
creasing monotonies lly . In practice the node numbers should be num- 
bered so as to minimize the bandwidth. Once the input is read the 
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program operates with the internal numbers which are now numbered 
from 1 through the number of nodes in the model. In this manner 
the node ordering and therefore the bandwidth of the stiffness 
matrix can be easily changed and nodes can be inserted or deleted 
by changing the external node specification.. 

The input is specified by entering the appropriate number on 
the input cards in fields of five. However, for this purpose the 
user can also utilize a shorthand form of the input. That is, spe- 
cifying m and -n consecutively is the equivalent of the speci- 
fication of nodes m, m 4- 1, m + 2, n and specifying m, -p, 

and -n consecutively is the equivalent of the specific tion of 
nodes m, m + p, m 4- 2p, ..., m 4- kp where m 4- kp is the highest 
integer of the form less than or equal to n. For example, the spe- 
cification of nodes 1 through 100 is written as 1 - 100 and 
1, 3, 5, ..., 99 as 1-2-99. This card input appears in fields of 
5 (15 Format) with lb items per card. Any number of continuation 
cards may be used. A zero or blank 15 field ends the specification . 

IV. Node Numbers for Members (A4,6X,4I5) 

The first field defines the element type: 

TRIR - triangular element designation (Ref. 1) 

SHEL - shell element designation (Refs. 2, 3) 

RING - ring element designation (Ref. 4) 

The first integer field designates a member number and the three 
following fields give node numbers for the member.. 

For triangular members, three nodes are specified. The i, j , k 
nodes are .given in counterclockwise order. For shell members, two 
nodes are specified, i and j, in the positive s direction. 

For ring members, only one node is specified. SEND in the first 
field ends the section . 

Note: Maximum number of elements (members) allowed is 600 currently 
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V . Global Coordinates (K15 . 7 , 1315) 

The coordinates are read in a specified order, r first, then 
z, and last the PHI coordinates. PHI is the angle in radians between 
the z-axis and the normal to the shell surface, measured positive 
clockwise (Fig. 2) from the z-axis. The floating point number desig- 
nates a coordinate value. Following are the applicable node numbers 
as in Section III. 

Only the integer fields are used on continuation cards. A zero 
or blank node field other than the first field on a noncontinuation 
card ends the card scan. On a continuation card, the floating field 
is left blank. A zero or blank first node field on a nonconti-nuation 

card ends each section for r, z, and PHI. If there are no shell 
elements, PHI is set to zero by default and need not be specified. 

V!. Boundary Conditions (611, 4X, 1415) 

The first six fields give a boundary condition specification, 
in the order: 

1) u r > u z > u q j if the node is attached to a triangular 
member. The last three fields for such a node must 
be zero. 

2 ) x, e g , v, Y g0 , if the node is attached to a 
shell element. See Fig. 2 for convention, u is the 
meridional displacement of the shell middle surface, 
w is the normal displacement, x is the rotation of 
a normal to the middle surface, e is the meridional 

b 

membrane strain, v is the circumferential displacement, 
and 7 g0 is the shear strain. At a pole, the condition 
u r -0, x-0, v = 0 , Tgg-0, must be prescribed. 
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3) u, vi, x if the node is attached to a ring member. 

The last three fields must be zero. 

Zero - denotes a fixed degree of freedom 

One - denotes a free degree of freedom 

Two - denotes an applied displacement and a corresponding 
card must be included in Section VII, designating 
the magnitude of the applied displacemen t . 

The integer fields are the applicable node numbers as in 
Section III. On continuation cards, only the integer fields are 
used. If a node’s boundary conditions are not specified in this 
section, all six degrees of freedom are assumed free. To change 
the default, the first card should be set to the desired default 
(i.e., for triangles 111000) with all nodes used in the problem 
then specified. 

Blank card ends section . 

Note: Maximum number of nodes currently allowed is 900. 

VII . Dependent Degrees of Freedom (2I5,2 (5X,2I5 < E15.7)) 

The first two fields designate a node number and its degree 
of freedom-. The dependency is defined in the following three 
fields. The two integer fields are node number and degree of 
freedom. The coefficient is contained in the following field. 

If there is another dependency for the node, it would be desig- 
nated in the last three fields . 

If there is yet another dependency, it would appear on a con- 
tinuation card with the first two integer fields blank. Both 

single and multipoint constraints of the form 6 i = a x 5 jl + 

R , -4- a 6 mav be considered as well as applied dis- 

2 j2 n jn 

placements of the form 5^ ~ a -y> e -S*> 

Blank card ends section. 
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Linking Up a Shell and a Triangle for a Cylinder 



Z> U i 


From shell theory restrictions 


u = n 

r 27 r 28 



w 


26 


u 

z 27 


u 26 ~ h/2 X 26 


u 

z 28 


u 


26 



U 26 + h/?2 X 26 
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The input cards, for this case, would appear as (let h = 0.1 inches) 


27 

1 

26 

2 

1.0 

E00 





28 

1 

26 

2 

1.0 

E00 





29 

1 

26 

2 

1.0 

E00 





27 

2 

26 

1 

-1.0 

E00 

26 

3 

- 0.05 

E00 

28 

2 

26 

1 

-1.0 

E00 





29 

2 

26 

1 

-1.0 

E00 

26 

3 

+ 0.05 

E00 


Note: In the boundary condition cards for nodes 27, 28, and 29, de- 

grees of freedom 1 and 2 must be fixed (specified = 0) 
because they are effectively eliminated from the solution 
(they are dependent degrees of freedom) . 

2 . Rotation of Displacements at Boundary 



For this problem, we set the displacement normal to the boundary 
equal to zero. We have: 

u, = u cos cp + u sin cp 

1 r z 

u„ = - u sin <p + u cos cp 

2 r z 
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Let us assume that the mode in question is node 25 and = 30 de- 
grees. The multipoint constraint cards are: 


25 

1 

25 

1 

0.86602540 E00 

25 

2 

0.50 

E00 

25 

2 

25 

1 

- 0.50 E00 

25 

2 

0.86602540 

E00 


Note: On the boundary condition cards, since u^ - 0, node 25 

degree of freedom 1 is specified to be 0 but node 25 
degree of freedom 2 is specified 1 (or fre^) since it 
remains an independent degree of freedom (although rotated) . 
If the normal displacement u^ were not fixed (and the 
transformation performed merely to apply a normaf load) it 
would have a 1 boundary condition specified. 

\ 

3. Applied Displacement 

Suppose for a shell we wish to apply a unit lateral (w) dis- 
placement at node 25. We have the multipoint constraint card as 

25 2 25 2 1.0 E00 

Note: A 2 -boundary condition for node 25 degree of freedom 2 

is specified on the boundary condition card. The end of 
this section is designated by a blank or zero in the third 
integer field (e.g., blank card). 
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VIII. 


Material Properties 


The format of this input is dependent upon the member type: 


A. 

Triangle 

Elements (4(5E15 .7) ,3E15 

7,/, (1615)) 


El 

Young's modulus in 1 direction 
(if £3 = 0, O-direction) 

Card 

E2 

Young's modulus in 2 direction 
(if (3=0, z-direction)' 

1 

E3 

Young's modulus in 3 direction 
(if £3 = 0, r-direction) 


V12 

Poisson's ratio in 

1-3 directions 


V23- 

Poisson's ratio in 

2-3 directions 


V31 

Poisson's ratio in 

3-1 directions 


G12 

Shear modulus in 

1-2 plane 

C'ard 

G23 

Shear modulus in 

2-3 plane 

• 2 

G13 

Shear modulus in 

3-1 plane 


S3 Orientation of material properties axis with 

respect to 0-z coordinates (in degrees) . 
This is a rotation about the r-axis and 
(3 is positive clockwise. 
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Card 

3 


Card 

4 


Card 

5 


Note : 


SIGOX Tensile (or compressive) yield stress 

in 0 -direction 

SIGOY Tensile (or compressive) yield stress 

in z-direction 

SXGOZ Tensile (or compressive) yield stress 

in r-direction 

SIGXY Shear yield stress in r-0 plane 

RMOSN Ramberg-Osgood shape parameter (n) , 

e = al E + (3o/7E) (ct/oq ^) n If set = 0.0 

material for the group of members is assumed to 
elastic ideally plastic. If RMOSS =0.0 and 
RMOSN # 0, RMOSN is ratio of slope of stress- 
strain curve for linear strain hardening to 
Young' s modulus (a = E^/E where Ey is the 
tangent modulus). 

RMOSS Ramberg-Osgood parameter cq 7 • If zero and 

RMOSN not zero have linear strain hardening. 

YLDST Yield' stress in tension (or compression) 

RMOSE Ramberg-Osgood parameter E (Young's modulus) 

Shear yield stress in 9-z plane 
Shear yield stress in z-r plane 


3 thermal coefficients of expansion in 
1, 2, 3 directions ( 0 , z, r if (3 = 0°) 


SIGYZ 

SIGZX 


TALPH (1) 
TALPH (2) 
TALPH (3) 


Orthotropic material properties should not be used for 

/ 

strain hardening problems in the current version of REVBY. 
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Succeeding cards give applicable members as in Section III. 

A zero or blank El ends the section. If there are no triangular 
members , no card is required . 


-B. Shell Elements ((4E15.7 ,/,5E15.7,/,3E15.7,/,2El5.7)/) , (1615)) 


1) Material Properties 


El 


Card 

1 


E2 

G33 


V12 


Young's modulus in meridional direction 
Young's modulus in circumferential direction 

Shear Modulus in plane of shell median surface 

Poisson's ratio in plane of shell median surface 


Card 

2 


SIGOX 

SIGOY 

SIGOZ 

SIGXY 

RMOSN 


Tensile or compressive yield stress in 
meridional direction 

Tensile or compressive yield stress in 
circumferential direction 

Tensile or compressive yield stress in 
direction normal to shell surface 

Shear yield stress in plane of shell median surface 
(See Card 3 for triangular members) 


RMOSS 
Card YLDST 

’ RMOSE 


(See Card 4 for triangular members) 


Card 

4 


TALPH (1) 
TALPH (2) 


Two thermal coefficients of expansion in meridional 
and circumferential direction, respectively 
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Succeeding cards give applicable members as in Section III. 

A zero or blank El ends the section. If there are no shell ele- 
ments, no blank card is required. 


2) Thickness for Shell Elements (follows the material 
properties. All nodes must have a thickness specified.) 

(E15 .7,/ (1615)) 

First field designates thickness value. Applicable nodes 
follow as in Section III. Section ends with zero or 
blank in first field. 

3) NG and NLRS (215). For shell elements, one card. NG 
designates the number of Gauss points used in integra- 
tion scheme for shell element, and the second field, 

NLRS, designates number of layers through the thickness 
at which plastic strains, stresses, etc., are calculated. 
NG < 8 , NLRS must be even, and < 30. 


C. Ring Elements (5E15 . 7 , / , 2E15 , 7 , / , A4. 6X, 3E15 . 7 , / ,3E15 . 7 , / , 
‘ 3E15. 7, /,E15. 7,(1615)) 


Card 

1 


E 

A 

IX 

IXY 


IY 


Young' s modulus 

Cross sectional area of ring 

Moment of inertia of ring about centroidal 
x'-axis (see diagram below) 

Product of inertia of ring about centroidal 
axis 

Moment of inertia of ring about centroidal 
y -axis 
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zo 


Card 

2 

ECCEN(Y) 


Distance from shell middle surface to ring 
centroid .(normal to shell middle surface) . 
Plus in negative z-direction . 

Distance along meridian from ring centroid to 
point of attachment to shell middle surface. 
Plus in negative s -direction. 


Card 

3 


SREC 
SC IR 
ZSEC 
I SEC 
HClR 


Solid rectangular cross section 
Solid circular cross section 
Z-section 
I-section 

Hollow circular cross section 


Remaining three fields designate section dimensions as shown in the 
diagrams below . 
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INPUT 


SREC 


Card T 2a 
3 < 2b 


SC IR 

ZSEC 

ISEC 

HCIR 

a 

a 

2a 

a 


2h 

2h 

t 


b 

b 



Card 

4 


{ 


t 


1 


t 


2 


t 


3 


Note: For ZSEC, a = b, t-^ = t^; for I SE C> 2a = 2b,' t^ = t^ • 


Card 

4 


If section name is ZSEC or ISEC, this card is required 
for three additional dimensions, otherwise this card 
is omitted. 


Card 

5 


SIGOX 
RMOSN ‘ 
KMOSS - 


Yield stress in tension or compression 
Same as shell and/or triangle 


Card TALPH (1) Thermal coefficient of expansion 

6 


Succeeding cards give applicable members as in Section III. 

A zero or blank E ends the section. If there are no ring elements, 
no card is required. 
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IX . Load Vector Components 


This format is also dependent on member type . 

A. Triangular Members (I5,3E15.7,/,I10,3E15.7) 

The integer field gives node number (i-j pairs); the floating- 
point fields give components at the respective nodes as follows: 

p =1.0 corresponds to an applied unit pressure 

(force/unit area) in the radial direction 
at node i 

p =1.0 corresponds to an applied unit pressure 
(force/unit area) in the z-direction 
at node i 

p e = 1.0 corresponds to an applied unit pressure 
(force/unit area) in the 0 -direction 
at node i 

The load between any two adjacent nodes is assumed to vary 
linearly. Therefore, the load vector is specified by two cards, 
the first specifying the load at node i and the second, the load 
at node j. A ze ro-or blank first integer field ends the section. 

If there are no uniform loads, a blank card is required. If there 
are no triangular elements , no card is required . 

B. Shell Members (3E15.7 ,/, (1615)) 

The three fields designate the meridional component, p , 

s 

normal component, p z , and the circumferential component, p^ , 
of applied stresses at nodes. Applied stresses are assumed to vary 
linearly from node to node. The applicable nodes follow as in 
Section III. A blank card ends the sec ti on. If the re a re no uni- 
form loads, a blank card is required. If there are no shell ele- 
ments, no card i s required. 
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C. Concentrated Loads (4E15 . 7 , / , (1615) ) 

For triangle component, specify F F , and F in r, z, 0 

1 Z L 

directions, and for the shell component specify N , Q , M , and 

o o o 

N in u, w, Xj and v directions, respectively, in lb/ in . 

These are forces per unit length of circumference, except at r - 0, 
where F = F = 0 and N = M = K , = 0 and F or Q repre- 

IT w S S S'7 Z b 

sent an actual force (lbs) . 

Applicable nodes follow as in Section III. A blank or zero 
card ends the concentrated loads specification section . 

X. Nodal Temperatures (A4,E15 . 7 , / (A4 ,3E15 . 7 , / , (1615) ) 

If the first field is SEND, there are no temperatures and this 
section ends. If the first field is blank, there are thermal loads 
and the floating number field designates the reference temperature 
value (T ref ) , which is only for the user's information. 

The next cards contain AT - T - T re £ at each node (up to 
three per node depending on the element to which it is attached) 
in the floating point fields. For a triangular element, one tem- 
perature is specified per node. For a shell element, three tempera- 
tures are specified, one at the top surface (+ h/2), T^, one at 
the median surface (z = 0), T 2 , and one at the bottom surface 
(- h/2), T^. The program assumes a parabolic distribution of tem- 
perature through the thickness (if T 2 = (T^ + T^)/2, this reduces 
to a linear variation). For a ring element, three temperatures are 
specified. A parabolic distribution is assumed through the ring 

rxi 

depth. T, is the temperature at x = +a or +h . T ? is the tern- 
perature at the centroid, and T^ is the temperature at x = -a or -h 
(see diagram in Section VIII) . 

Applicable nodes follow as in Section III. 
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XI. 


Members Whose: Rcsu3 L..- - ire to be Printed for Elastic Sol ufi.cn 
(163 5) /.s in Sect: ion I i I 

Blank cards end section. 

XII. Nodes Whose Results n^e to be Printed for Elastic Solution 
I16I5L As in Section III 

Blank cards end section. 

XU. Same as Section XI for Plasticity, if applicable. 

Blank cards end section. 

XIV. Same as Section XII for Plasticity, if applicable. 

Blank cards end section. 

XV . PMAX and PPCT (2E15-7) 

PMAX - Maximum load to be applied for this half cycle 
PPCT - Load increment as a percentage of yield load 

XVL NPRNTj, PMAX, PPCT (I5.2E15.7) 

For succeeding load cycles, one card giving now NPRNT, PMAX, 
and PPCT. Zero NPRNT signifies no new load cycle and end of prob- 
lem. For succeeding half-cycles, PPCT equal to one-half the original 
value should bemused since the elastic range is twice the elastic 
range for initial loading. 

XVH. Change in Material Properties (A4) 

To change properties of any group of members, read indicator 
for new material properties (A4) . 

A. If RING is input, it indicates a change in the RING element 
material properties. A card with a member number belonging 
to the group whose properties arc to be changed and the new 
quantities arc read in (I5,4E15 . 7) . 
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Card 

1 


Card 

2 


Card 

3 


C . 


Member Number 

SIGOX 

RMOSN 

RMOSS 

TALPH 


Member of group to be changed 


New Values 


NEWM signifies new material properties for either triangle 
or shell elements. The new values are read in (I5,5E15.7 ,/, 
5E15 . 7 j / , 3E15 . 7) 

Member Number Member of group to be changed 

SIGOX 

SIGOX 

SIGOZ New Values 

SIGXY 

RMOSN 

RMOSS 
YLDST 

RMOSE New Values 

SIGYZ 

SIGZX 


TALPH 

(1) 


TALPH 

(2) 

New Values 

TALPH 

(3) 



SEND as input indicates no new material properties or end of 
this section. 
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XVIII. Each Problem* s Input Must Be Ended with a Card Reading 

"ENDb" Where "b" Denotes a Blanks in Columns 1-4. The 
Last Problem in a Run Should End with a Card Reading 
"STOP" in Column 1-4 Instead of "ENDb." 
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SECTION 5 


INSTRUCTIONS FOR USE OF 
OUT -OF -PLANE GM 

A Program for the Nonlinear 
Analysis of Built-Up Structures 




OUT-OF-PLANE GM is a finite element program for the nonlinear 
analysis of arbitrary built-up thin walled structures. This pro- 
gram considers the sheet material to be' in states of membrane 
stress. The finite element library consists of the following 
elements : 

9 Three-node uniform stress triangle 

® Six-node linearly varying stress triangle 

® Four- and five-node hybrid triangles to be 
used as transition elements between three - 
and six-node elements (see Fig. la) 

6 Two-node uniform stress stringer 

* Three-node linearly varying stress stringer 

® Beam with various cross sections subjected to 
bending about two planes as well as torsion 

The program _is_ capable of treating the elastic and the elastic- 
ideally plastic response of orthotropic materials.. In addition, con- 
sideration is given to isotropic materials exhibiting elastic-ideal. 
plastic, linear strain hardening, or nonlinear strain hardening be- 
havior. Further, the kinematic hardening theory of plasticity i. used. 

The effects of geometric nonlinear behavior is treated by c<. n- 
sidering a converted coordinate system that accounts for changes in 
geometry and by the introduction of an initial stress stiffness matrix. 
The solution strategy implemented by the program is outlined in a com- 
panion volume, "PLANS - A Finite Element Program for Nonlinear Analy- 
sis of Structures, Volume I - Theoretical Manual," NASA Contractor 1 
Report NASA CR-2568, p. 31. 

The input to the program is categorized in the following sections. 
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Any 80-character title describing the problem. 



CARD 1: 

NPNTC is the sum of the following integers corresponding to the 
option desired. 

If NPNTC =0 No intermediate printout 

= 1 Print the load vector 

= 2 Print element stiffness matrix 

= 16 Print each element stiffness matrix entry to be 
stacked with its stacking index 

= 32 Print the total stiffness matrix 

For example, if it is desired to print the load vector and the 
total stiffness matrix, NPNTC = 1 + 32 = 33. 

NPRNT > 0 Print output every NPRNT increments of load. 

If NPRNT equals zero the output is printed 
every load step. 

NFORM 2 ' 0 Reform elastic stiffness matrix to account for 

changing geometry every NFORM increments of load. 
The effect of the geometric (initial stress) 
stiffness matrix is still introduced as an "ef- 
fective load" in this case. 

. IREST: IREST =0 Do not generate a restart tape 

IREST = 1 Build a restart tape 
IREST = 2 Run continuing from a created 
restart tape 
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NUTAP: 


INPRT : 


NRSRT : 


CARD 2: 

PMAX: 


DELP: 

PTAN: 


Applicable only if restarting from a restart 
tape. When IREST = 2 and NUTAP = 0 no additional 
restart data is written. If NUTAP ■ 0 restart 
tape is continued. 

INPRT = 0 Write restart tape only at P - PMAX 
(i.e. , at maximum load) 

INPRT = N Restart data written every INPRT 
time steps 

Applicable only if restarting from a restart tape. 
NRSRT is a unique number that is printed in the 
job that generated the restart tape for each load 
that data is written on the restart, tape. This 
number defines the starting load for the continua- 
tion job (see Section XVII). 


Maximum load factor. The load factor multiplies 
the vector of applied load. If the applied load- 
ing is normalized to unity then PMAX represents 
the magnitude of the largest component of the 
plied load. If the maximum desired load is input 
then PMAX is equal to one and each load step repre 
sents a percentage of the total load. 

Number of load increments. The load increment 
AP is specified by PMAX/DELP. 

Load at which full tangent modulus method is used, 
i.e., the stiffness matrix explicitly includes the 
initial stress stiffness matrix. If the "tangent 
modulus" method is to be used for the entire a'>\- 
sis, PTAN = 0.0. In this case the stiffness P’.-tiix 
is reformed in every load step. If the "effective 
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load" method is to be used for the entire analy- 
sis, PTAN - PMAX. In this case the stiffness 
matrix includes only the elastic stiffness matrix 
and is reformed every NFORM load steps. When 
0 < PTAN < PMAX then the program uses the "ef- 
fective load" method up to PTAN and then imple- 
ments the "tangent modulus" until the maximum 
load is reached. 

EQMULT: Factor multiplying the equilibrium correction 

term. EQMULT =0.0 indicates no equilibrium 
correction and EQMULT =1.0 indicates full 
equilibrium correction. 

III. Node Specification (1615) 

This section defines an allowable set of external node point 
numbers. The maximum node number that can be used is 9999. The 
program uses this information in two ways. First to set up a table 
of allowable node points that is used to check all subsequent node 
point input. Secondly, the program converts each external node num- 
ber to an internal number consecutively in the order that the node 
appears on the input card. Consequently the order of the input of 
external node numbers is completely arbitrary and need not be in- 
creasing monotonically . In practice the node numbers should be num- 
bered so as to minimize the bandwidth. Once the input is read the 
program operates with the internal numbers which are now numbered 
from 1 through the number of nodes in the model. In this manner 
the node ordering and therefore the bandwidth of the stiffness 
matrix can be easily changed and nodes can be inserted or deleted 
by changing the external node specification. 

The input is specified by entering the appropriate number on 
the input cards in fields of five. However, for this purpose the 
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user can also utilize a shorthand form of the input. That is, spe- 
cifying m and -n consecutively is the equivalent of the speci- 
fication of nodes m, m + 1, m -h 2, ...,n and specifying m, -p, 
and -n consecutively is the equivalent of the speeffiction of 
nodes m, m + p, m + 2p, m + kp where m + kp is the highest 

integer of the form less than or equal to n. For example, the spe 
cification of nodes 1 through 100 is written as 1 - 100 and 
1, 3, 5, 99 as 1-2-99, This card input appears in fields of 

5 (15 Format) with 16 items per card. Any number of continuation 

cards may be used. A zero or blank 15 field ends the specification 

IV . Member Connectivity (Node Number of Each Member) 

FORMAT (A4, fiX, 915) 

The first alphanumeric field defines the element type; 

TRIM - Triangular membrane element (Ref. 3) 

■ BEAM - ! Beam element 

STRG - Stringer element 

The firsf integer field designates the member number. The next 
integer fields designate the connecting nodes as follows: 

TRIM - The nodes for the triangular family of elements 
are specified around the perimeter beginning with a maior (vertex) 
node. Then a minor (midside) node and then alternatively maior and 
minor as shown in Fig. 1. The absence of a minor node must be in- 
dicated by a zero or blank ^ield in the proper position. 

BEAM - Three node specifications are necessary- for the 
beam element. 
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Uniform Strain 
Elements 




Pig. la Triangular Family of Finite Elements Used 



Fig. lb Elements Topology 
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Nodes i and i (Fig. 2) which designate the element end 
points and a third node k, defining the normal to the beam axis 
about which the section properties are specified. This additional 
node may be a node of the structural idealization or it may be a 
"fictitious node" specified lust for the purpose of defining the 
beam section properties. This is shown in Fig. 2. 

STRG - Three node specifications are necessary for the 
stringer element. Nodes i and j, connecting the end points, 
and if desired an additional node designating a midpoint node. 

This is shown in Fig. 3. A zero or blank for the midpoint node 
specifies a two-node stringer. The midside node is the third node 
specified. 

SEND Ends the section 

V. X -Coordinates of Nodes FORMAT (E15.7, 1315) 

The x-coordinates of the nodes appearing in the 15 fields 
are set to the value in the E15.7 field. Any number of continua- 
tion cards may be used; their first fifteen columns are ignored. A 

zero or blank 1-5 — field terminates the card scan for a given x- 
coordinate. A zero or blank first 15 field (columns 16-20)' on a 
noncontinuation card ends the section. Both shorthand representa- 
tions of Section III are allowed. 

VI. Y-Coordinates of the Nodes . Same as Section V. 

VII. Z -Coordinates of the Nodes . Same as Section V. 

VIII. Boundary Conditions FORMAT (611, 9X , 1315) 

The first six fields give the boundary conditions specifica- 
tions in the order: u, v, w, a , P , P ; where u, v, w are 

x y z 

global displacements in the x, y, z directions, respectively, 

and 0 , P , 9 are rotations with respect to the x, y, z axes, 
x y z J 

The last three fields for , P , 9 are only used when using a 

x y z 

beam element. 
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Zero 


denotes a fixed degree of freedom 
One denotes a free degree of freedom 

Two will result in the application of a unit general- 
ized displacement, or a corresponding card may be 
included in Section IX (dependent degrees of free- 
dom) designating the magnitude of the applied gen- 
eralized displacement. 

The 1315 fields give the applicable nodes for the boundary 
condition specifications, with both shorthand notations of Sec- 
tion III permitted. Any number of continuatuion cards may be used 
for a given specification. However, only the 1315 fields are 
used on a continuation card. A zero or blank 15 field terminates 
the card scan for a given boundary condition specification. Note : 

If the last field of a card (columns 75-80) is the last specifica- 
tion, an additional blank card (continuation card) must follow. A 
zero or blank first 15 field (columns 15-20) on a noncontinuation 
card ends the section. If a node’s boundary conditions are not 
specified in this section, all the degrees of freedom are assumed 
to be free. To change this default condition, the first card of 
this section should be set to the desired default (i.e., in the 
absence of beams, 111000) with all nodes used- in the problem spe- 
cified. Note: Maximum number of nodes is currently 900. 

IX . Dependent Degrees of Freedom FORMAT (215, 2(5X, 215, E15.7)) 

This section designates the input for both single and multipoint 
constraints as well as applied displacement of the form: 


1) 6. = a.6... + 5 . 0 + ... a 6 

' x 1 jl 2 i2 n n 


where 5^ is a dependent degree of freedom, 

5^ ... are independent degrees of freedom 

and a, , . . . , a are coefficients 

12 n 
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2) Rotation of displacements at a node 



A/ 

A> 


A t 

6 i 

= a, 6 . 

1 i 

4- a„S . 
2 1 

+ 

a 3 5 k 


Aj 

A* 


A> 

cn 

L_u 

“ ^l 6 i 

+ (3 2 6. 

+ 

P 3 6 k 


a) 

Aj 


A» 

5 k 

= ■y 1 8 1 

+ 


Y 3 6 k 


•where the a, P, v’s are the direction cosines of the 
- rotation, 5., 6., fi, are the displacements with re- 

1 1 K /\j s\_, 

spect to the original global directions, and 5., 5., 6, 

1 J K 

are the components of displacements at the node with 
respect to the new coordinate axes. An example of this 
capability is given in Appendix I. 

3) Applied generalized displacement 



where the coefficient is the applied generalized 

displacement . 

The first two fields designate a node number and a degree of 
freedom (i.e., 1-6). The dependency is defined in the following 

three fields. The two integer fields designate the node number 
and degree of freedom number and the coefficient by the floating 
point field. If there is another dependency for the node, it is 
designated in a similar fashion in the next three fields. Any 
number of continuation cards can be used with the first two fields 
blank. The section is ended by a blank or zero in the third in- 
teger field (blank card) . Examples of the use of multipoint con- 
straints are in Appendix I. 

X. Material and Section Properties 

The format for this input is dependent upon the member type. 
Each type of input begins with a word of up to four characters. 
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MATl Material properties for plane stress, necessary 
with membrane triangles (TRIM) 

MBM Mate rial properties for a beam element 

MSTG Material and section properties for a stringer 

THIK Member thickness for triangular membrane (TRIM) 

* 

elements 

MBET Angle between local axes and principal directions 
of orthotropy for TRIM elements 

SREC Beam section properties for a solid rectangular 
section 

SCIR Beam section properties for a solid circular 
section 

ZSEC Beam section properties for a Z-section 

ISEC Beam section properties for an I-section 

HCIR Beam section properties, for a hollow circular 

section 

HREC Beam section properties for a hollow rectangular 
section 

LSEC Beam section properties for an L-section 

TSEC Beam sectiqn properties for a T-section 

CSEC Beam sectiqn properties for a channel section 

SEND Ends the section 

MATl - Plane Stress Material Properties - FORMAT (A4, IX, 5E15.7, /, 

2E15.7, /, 4E15.7, /, 
5E15.7, /, (1615)) 

The first four cards specify material properties, as follows: 
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CARD 1: MAT1 


EONE = Young's modulus in principal property axis (1) 
ETWO = Young's modulus in principal property axis (2) 
BETA = No longer used. Set equal to zero. 

GONTO = Shear modulus in (1) — (2) principal property 
plane 

VONTO = Poisson's ratio, v ^ 


Note: 





y 12 



CARD 2; TALF-1 - Coefficient of thermal expansion in 

1- axis direction, a ^ 

TALF-2 = C oefficient of thermal expansion in 

2 - axis direction, 

CARD 3: SIGOX = Yield stress in principal l~direction 

SIGOY = Yield stress in principal 2-direction 

S1G0Z = Yield stress in principal 3-direction 

SIGXY = Shear yield stress in principal 1-2 plane 

CARD 4: RMOSN = If RMOSS ^ 0; RMOSN = n, the shape parameter 

used in Ramberg-Osgood representation of 
stress-strain behavior 


If RMOSS = 0; RMOSN = a the slope of the 
linear strain hardening stress-strain repre- 
sentation, i.e., a = E^/E where E^ is 
the tangent modulus 
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RMOSS = If RMOSS ^ 0; RMOSS = Ramberg- Osgood 
parameter <Jq ^ 

Note 1: If RMOSN = 0 and RMOSS - 0, 

the material for the element(s) 
is assumed to be elastic-ideally 
plastic 

RMOSE = Ramberg-Osgood parameter E (Young’s modulus) 

YLDST = Yield stress in tension 

YLDSC = Yield stress in compression 

Input for the yield stress in tension and compression has been main- 
tained in order to accommodate materials that exhibit initial aniso- 
tropic plastic behavior. In this case an initial translation of the 

u. 

1 

yield surface is made consistent with the kinematic hardening theory. 

Note: Only initially isotropic materials can be treated when con- 

sidering linear or nonlinear strain hardening using this 
property card. 

Succeeding cards- give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing . 


MBET ~ Orientation of Axes of Material Anisotropy - FORMAT (A4, IX, 

E15.7, (, (1615)) 


CARD 1 : MBET 

BETF = Angle ('• in degrees between local x-axis and 
principal 1-axis for material orthotropy. See 
Fig. 4. Only applicable for TRIM elements. 


l 

Isakson, G. , Armen, H. Jr., and Pifko, A., "Discrete-Element Methods 
for the Plastic Analysis of Structures," NASA Contractor's Report NASA 
CR-803, October 1967. 


5-13 



Note: 


This card is an optional card. The de- 
fault is BETF =0.0 for all members. 
It should be used with initially ortho- 
tropic materials. 


x 


Fig. 4 Orientation of Element Local and Material Axes 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero ar blank 15 field ends 
each member listing. 

MSTG - Stringer Properties - FORMAT (A4, IX, 5E15.7, /E15.7/, (1^15)) 

CARD 1 : MSTG 

E = Young’s modulus 

A - Cross sectipnal area 

RMOSN = If RMOSS 4 0; RMOSN = n, the shape parameter 
used in Ramberg-Osgood representation of 
stress-strain behavior. 
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If RMOSS = 0; RMOSN = a, the slope of the 
linear strain hardening stress-strain repre- 
sentation, i.e., a = E^/E where E^ is the 
tangent modulus 

RMOSS = If RMOSS ^ 0; RMOSS = Ramberg-Osgood 
parameter cTq j 

Note 1: If RMOSN = 0 and RMOSS =0, the 

material for the element(s) is 
assumed to be elas tic-ideally plastic 

YLDST = Yield stress 

CARD 2; TALF = Coefficient of thermal expansion 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted- Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing- 


THIK - Element Thicknesses - FORMAT (A4, IX, E15.7) 

Necessary with membrane elements 
TRIM 


CARD 1 : THIK 

THICK = Element thickness 


Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing; . 


MBM - Beam Material Properties - FORMAT (A4, IX, 5E15.7, /, Elf. 7, 

(1415) ) 


CARD 1 : MBM 

E = Young 1 s modulus 
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ANU = Poisson 1 s ratio 

RMOSS = If not equal to zero, RMOSS equals 
Ramberg- Osgood parameter, cr^ ^ 

RMOSN = If RMOSS ^ 0; RMOSN - n, the shape parameter 
used in Ramberg- Osgood representation of stress- 
strain behavior. 


If RMOSS = 0; RMOSN = a, the slope of the 
linear strain hardening stress-strain repre- 
sentation, i.e., a = E t /E where E^ is the 
tangent modulus 

YT.DST = Yield stress 

CARD 2 : TALF = Coefficient of thermal expansion 


Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing . 

Beam Section Properties - FORMAT (A4, IX, 5E15.7, /, 3E15.7, /, 

formats for Cards 3 and 4 (see 
below), /, (1615)) 

The cards specifying beam section properties all start with 
the following information: 


CARD 1: 


SRECh 

SCIR 

ISEC 

ZSEC 

HCIR \ 

HREC 

LSEC 

TSEC, 

CSECy 


Defined at the beginning of the section. 
The different cross sections are shown 
in Tables 1 and 2. 
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A = cross-sectional area 

Iyy = moment of inertia about y axis (see Fig. 5) 

Izz = moment of inertia about z axis (see Fig. 5) 

Iyz = product of inertia 

J = torsional rigidity 

CARD 2: YO = Eccentricity of attachment point in the y^ 

direction 

ZO ~ Eccentricity of attachment point in the 
direction 

[3 = angle defining the transformation of the 

ly, Iz, Iyz to another reference axis 
(see Fig. 5) 

Additional cards are required according to which section is speci- 
fied. The notation for each section is shown in Table 1. 



Fig. 5 Definition of Coordinate Axes in 
Cross Section of Beam Element 
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SREC - FORMAT (2E15.7) 

CARD 3: A = Width 

B = Depth 

SCIR - FORMAT (E15.7) 


CARD 3: 

R 

= Radius 



ZSEC, ISEC, 

and 

CSEC - FORMAT 

(3E15.7, 

/ , 3E15.7) 

CARD 3: 

A1 

= Dimension 

of upper 

flange 


A2 

= Dimension 

of web 



A3 

= Dimension 

of lower 

flange 

CARD 4: 

T1 

- Thickness 

of upper 

flange 


T2 

= Thickness 

of web 



T3 

= Thickness 

of lower 

flange 

HCIR - FORMAT (2E15.7) 

* i 



CARD. 3: 

R 

- Outer radius . 



T 

* = Thickness 



HREC, LSEC, 

and 

TSEC - FORMAT 

(4E15.7) 


CARD 3: 

Al 

= Width (parallel to 

y-axis) 


A2 = Depth (parallel to z-axis) 

T1 = Thickness of upper and/or lower flanges 
T2 = Thickness of vertical webs 

Succeeding cards give applicable members; both shorthands of Sec- 
tion III are permitted. Any number of continuation cards may be 
used for a given specification. A zero or blank 15 field ends 
each member listing. 

SEND Ends the section. 
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TABLE 1 SECTION INPUT AND GEOMETRY 




os-g 


TABLE 1 SECTION INPUT AND GEOMETRY (C.ONT) 


c . 

Section 

Keyword 

Input 

Hollow 

Rectangular 

Section 

HREC 

Area, moment of inertia, I , I , I , 

’ yy' zz yz’ 

J eccentricities, y , z , transformation 
angle p. Width, a^, depth, a^, thickness of 
upper and lower flanges, t^, thickness of 
vertical webs, 

L- Sect ion 

LSEC 

Area, moment of inertia, I , I , I , 

’ yy zz’ yz’ 

J eccentricities, y , z q , transformation 
angle p. Dimension of flange and web a^,a 
thickness flange, and web, t^, t^. 

T" Sect ion 

TSEC 

Same as for Ir Section 
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TABLE 1 SECTION INPUT AND GEOMETRY (CONT) 


Section } Keyword 


Input 


Z-Section ZSEC Area, moments of inertia, I , I , I , 

yy zz yz 

J eccentricities, y.,z, transformation 
angle . Dimension of upper and lower 
flange, a a^.) Dimension of web , > 

thickness of upper and lower flange, t^, t^, 
thickness of web, 


NOTE : 


a 3’ C 1 " £ 3 


T 1 \ 2 V 




il-Section | ISEC jSame as Z- section 

a l ^ a 3 ’ fc i ^ 


CSEC {Same as Z“Section 


a l a 3 ; t l C 3 



H~ A 3"H 
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TABLE 2 NUMBER AND LOCATION OF STRESS 


Keyword Number 
Stress 
Points 


Location of 
Stress Points 


Solid 

Circular 


± y 

0.430568a 
0. 430568a 
0.169991a 
0.169991a 


+ z 

0.430568b 
0.169991b 
0.430568b 
0. 169991b 



Location of 
Stress Points 



Shear Center at Centroid 



I Shear Center at Centroid 
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TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 


Keyword Number 
Stress 
Points 


Location of 
Stress Points 


352.8 


0. 77459a, 


0. 5a. 


0.5a. 


323.4 

36.6 

274.6 

85.4 

213.0 

147.0 


0.5a, 


0.5a, 


0. 77459a, 


Location of 
Stress Points 



Shear Center at Centroid 



a l * a l ‘ C 2 


Shear Center at Centroid 


a 2 a 2 " fc l 
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TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 
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TABLE 2 NUMBER 


Section 

Keyword 

Number 

Stress 

Points 


Z-Section 

ZSEC 

9 

1 

2 

3 

4 

5 . 

6 
7 

. 8 
9 

a 2 = 
a 2 


Location of 
Stress Points 


Location of 
Stress Points 


y 

-0. 88730a^ 
-jo. 11271a 1 
-O'.SOOOa^ 

0. 88730a 1 
0. 1127la 1 
0.5000a 1 * 
0.0 ' 

0.0 

0.0 

(a 2 - t 1 )/2 

( a 2 " V 


z 

a 2 
*2 
. *2 
*2 
a 2 
a 2 

0.77459a' - 

0.77459a' 

0.0 




Shear Center at 


. 

il 
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TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 


Section j Keyword 


CSEC 


Number 

Stress 

Points 


Location of 
Stress Points 


0. 88730a, -t- e 0.5a, 


0. 11271a, + e 0.5a, 


0. 5000a, + e 0. 5a, 


0. 88730a, + e ~0.5a. 


0.11271a, + e -0.5a. 


0.5000a. +e -0.5a, 


Location of 
Stress Points 



0. 77459a, 


Axes at Shear Center 


-0. 77459a, 


e - distance from shear center to 

g 

■web center line 


e = 
s 


l/2(a 1 - t 2 ) 


1 -F 1/6 


*2< a 2 ~ *!>] ’ a 2 =a 2" t l 
t l^ a l ‘ t l\ 


NOTE: t l = t 3 








\ 

N> 


TABLE 2 NUMBER AND LOCATION OF STRESS POINTS (CONT) 










XI. Applied Loads 


Two different types of loading are currently admissible 
concentrated loads at nodes, and distributed line loads on beams. 

Each loading situation is designated by one of the following 
four- character key words: 

CONC (Concentrated force or moment at a node) 

BMLO (Distributed line load, beam element) 

SEND (Section end) 

Note; The keywords for the section on applied loads are used 
as section headings. The keyword appearing on an 'input 
card designates that the input to follow is associated 
with a particular type of applied loading. A blank card 
(as described in succeeding sections) delimits the input 
and indicates that the next card contains a different key- 
word. For example, CONC indicates input that follows is 
for concentrated loads and TRIA; distributed loads on an 
edge of a membrane triangle. Thus the input is as follows: 

CONC 

data 

for 

concentrated 

loads 

blank card 

BMLO 

SEND 

It should be noted that the keywords can appear in any order 
and may be specified more than one time in the applied load 



CONC - Concentrated loads - FORMAT (A4, /, 15, 3E15.7, /, 5X, 3E15.7) 

Each card gives the load components in global directions at a 
specified node. 

CARD 1: CONC 

CARD 2: 15 field: Node number 

First E15.7 field: Force component 

Second E15.7 field: Force component F^ 

Third E15.7 field: Force component F^ 

CARD 3: First El 5. 7 field Moment component M^ 

Second E15.7 field Moment component M^ 

Third E15.7 field Moment component 

A blank card (i.e., zero or blank first 15 field) ends the speci- 
fication of concentrated loads. 

BMLO - Distributed Line Load, Beam Element - FORMAT (A4, /, 4E15.7) 
CARD 1 : BMLO 

Force/unit length in local y-direction (see 
'Fig. 5) at node i 

Force/unit len.gth in local y-direction at 
node j 

Force/unit length in local" z-direction at 
node i 

Force/unit length in local z-direction at 
node i 

Note: A linear variation of the distributed 

load between nodes is assumed. 
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Succeeding cards give applicable members; both shorthands of Section III 
are permitted. Any number of continuation cards may be used for a given 
specification. A zero or blank 15 field ends each member listing. 

A card with SEND in the first four columns ends the section for 
applied loading . 

XII. Members to be Printed — FORMAT (16I5 1 ) 

Specify the members whose strains and stresses are to be printed. 
Both shorthands of Section III are allowed. A maximum of 600 members 
may be specified. Members in excess of 600 and undefined member 
numbers are ignored. A blank card or card with only zero entries 
ends the section. 

XIII . Nodes to Be Printed — FORMAT (1615) 

Up to 900 nodes whose displacements are to be printed for the 
analysis, as per Section XII. 


XIV . Parameters for Succeeding Load Cycles — FORMAT (315, /, 
3E15.7) 

CARD 1: NPRNT = If equal to zero, no additional load cycle, end 

of problem. If greater than zero, print output 
every NPRNT increments of load 
NFORM = Reform elastic stiffness matrix every NFORM 
increments of load 

NNEWL = If equal to zero, no new load distribution is 
to be input. If greater than zero, a new load 
distribution will be read in 


CARD 2; PMAX - New maximum load factor 

DELP = Number of load increments for new load cycle 
PTAN = Load at which full tangent modulus is used 
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XV . New Applied Loads 

New load distribution is read in if NNKWb ^ 0. Formal is 
identical to Section XI. 

XVI. Problem End - FORMAT (A4) 

Each problem's input must be ended with a card reading "ENDb" 
where "b" denotes a blank, in columns 1-4. The last problem in 
a run should end with a card reading "STOP" in column 1-4 instead 
of "ENDb." 
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XVII. 



The Initial restart tape is created on Unit 21. Subsequent 
restart jobs mount the restart tape as Unit 21 and if desired 
.(NUTAP v 0) copies and continues the restart tape on Unit 22. 


CARD 1: Problem Title - FORMAT (2A4) 


As in Section I 


CARD 2: NPNTC, NPRNT , IFORM, IREST, NUTAP, INPRT , NRSRT 

As in Section II. Here NPNTC is ignored and 
IREST = 2 


NRSRT is a unique number obtained as output from the 
job that generated the restart tape. This output is 
of the form: 

,* * * * * 

RESTART TAPE HAS BEEN WRITTEN FOR P - 7 . 100000E+04 
NRSRT = 8 . 

***** 

In this case NRSRT is equal to 8 and the job will be 
restarted at a load factor of 71.0. A restart job 
can be started at any intermediate load level for 
which the restart tape has been written. 

CARD 3: PMAX, DELP, PTAN, EQMULT 

As in Section II 

Succeeding cards follow as in Sections XII through XVI. 
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APPENDIX I 


MULTIPOINT CONSTRAINT S 


EXAMPLE 1 



Node I and 2 are fixed. 

There is an applied displacement: at node 3 in the negative y-direction. 
The distance between nodes 3 and 4 remain constant. 


Boundary Conditions 

011000 1-4 default condition 

000000 1 2 4 

021000 3 

Blank card 


Dependencies 

Node Dof 

3 2 

4 2 

4 3 


Node Dof 

3 2 0.050 

3 2 1.0 

3 3 1.0 


(specified displacemei ' ) 
restraint condition that 
3-4 remain riftid 


NOTE: In the boundary condition cards for node 4, decrees of 

freedom 2 and 3 must be specified 0 because they are 
effectively eliminated from the solution (they are de- 
pendent decrees of freedom) . 
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EXAMPLE 2 


Rotation of Displacements at a- Boundary 


x 
25 

For this problem we set the displacement u equal to zero and 
v free. This corresponds to the normal displacement along x 
equal to zero and the tangential displacement (in y direction) 
free. The global coordinates of thepproblem are x, y. Thus: 

u = u cos 9 - v sin 9 

v = u sin 9 4- v cos 9 

The multipoint constraint cards are: 

25 1 25 1 0. 86602 540E00 25 2 -O.5OEO0 

25 2 25 1 0.50EOO 25 2 0. 85602 54E00 

Note: On the boundary condition cards, since u <= 0.0 node 2 5 

degree of freedom 1 is specified to be 0 but node 25 

degree of freedom 2 is specified 1 or free since it re- 

mains an independent degree of freedom (although rotated) . 
If the normal displacement u were not fixed (and the 
transformation performed merely to apply a normal load) , 
it would have a 1 boundary condition specified. 
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SECTION 6 


INSTRUCTIONS FOR' USE 
OF SATELLITE 


A Preprocessing Program 
for PLANS 



INSTRUCTIONS FOR USE OF SATELLITE 


SATELLITE is the preprocessing program associated with the 
PLANS system programs that checks and plots the undeformed finite 
element model. Additionally, it calculates the ordering of the 
external node numbers that leads to the minimum semibandwidth. 
SATELLITE is currently not operational for the HEX program. 

The input for SATELLITE is as follows: 

CARD 1: FORMAT (20A4) 

Keywords left justified in fields of four. One of the fol- 
lowing is required and indicates the program whose data follows * 

PLNE - PLANE program, In-Plane Membrane Program 
Used for Fracture Analysis 

BEND - BEND program, Elastic-Plastic Analysis of 
■Built-Up Structures 

REVB - REVBY program, Elastic-Plastic Analysis of 
Bodies of Revolution 

OPGM - OUT-OF-PLANE GM, Nonlinear Analysis of 
Built-Up Structures 

DYCA - DYCAST program, Nonlinear Dynamic Analysis 
of Built-Up Structures 

Additional optional keywords can be specified in fields of- four. 
Any number of spaces can be skipped on the input card that are 
multiples ol four. 

BAND - Specifies that the optimum semibandwidth is 
to be calculated. In this case the optimum 
order of the external node numbers for 
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Section III data will be printed and punched. 
These punched data are used with the input 
deck. If this keyword is not specified the 
bandwidth will not be optimized. 

SCAN - Process and check input data but do not 
produce any plots 

PRIN - Print transformed coordinates to be plotted 

The input data deck follows Card 1. Input through the nodal co- 
ordinates is necessary. The remainder of the deck may be input 
although the cards are read and not processed until a STOP or END 
card is reached. Input for the SATELLITE program continues as 
follows . 


CARD 2: FORMAT (515) 

This card determines the number and type of plots. Up to five 
different labeled pictures can be obtained. These are specified 
with the fixed point numbers right justified in any of the fields 
of five as follows: 

1 - labeled unconnected nodes 

2 - labeled nodes, unlabeled members 

3 - labeled nodes and members 

4 - unlabeled nodes, labeled members 

5 - un.labeled nodes and members 


CARD 3: 


FORMAT (4E15.7) 

ALPHA - rotation of structure about the global x-axis 

BETA - rotation of structure about the global y-axis 
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GAMMA - rotation of structure' about the global z-axis . 

LENGTH - maximum width of picture in the paper trans- 
verse direction. LENGTH is a function of the 
plotter being used. 

If ALPHA, BETA, GAMMA are zero then the picture is projected on the 
x-y plane with the y direction in the transverse paper direc- 
tion. The angles lead to a path dependent- transformation by ro- 
tating the body first abodt the global x then y and z axes, 
respectively. The resulting viewing plane is still the original 
global x-y plane. The viewing plane can be changed with the 
following card. 

CARD 4: FORMAT (6Al) 

This card .specifies the viewing plane to be used by specifying 
the six alphanumeric characters XYZ+-H-. The first two characters 
specifying the viewing plane with the first indicating the length 
direction and the second the width of the paper. The last three 
characters can be the characters + or - and indicate that the 
normal to the plane is in the positive or negative coordinate 
direction. For example: 

YZX+4+ indicates the viewing plane as the global YZ 
plane with the Z axis in the paper width 
direction and Y along the paper length 

CARD 5: FORMAT (4(A4,1X)) 

Element types to be omitted. For BEND, OPGM, DYCAST: 

BEAM - will omit all beam elements 

TRIM - will omit all triangle membrane elements 
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STGR will omit all stringer elements 

TRIP - will omit all triangular plate elements 

For REVBY: 

TRIR - will omit all triangle elements 

SHEL - will omit all shell elements 

RING - will omit all ring elements 

If there are no element types to be omitted a blank card must be 
specified. 

CARD 6: FORMAT (1615) 

Nodes to be included in the plot. All the shorthand notations 
are allowed, i.e., 1 through 100 is input as 1 - 100 and 
1, 3, ..., 99 as 1-2-99. A blank card or card with only zero 
entries ends the section . 

If only a part of the structure is to be plotted, this section 
should include only the nodes that lie in the section. Only these 
nodes are considered_wben scaling the plot so that the subsection 
will be scaled up to the maximum allowable size. 

CARD 7: FORMAT (1615) 

Elements to be included in the plot. Input is the same as 
for Card 6. All nodes associated with the elements to be plotted 
must be specified with Card 6. 

Note: If all the elements and nodes are to be plotted it is 

sufficient to specify one negative number in the first 
15 field whose magnitude is greater than the largest 
external element/node number in the finite element grid. 


6-4 



CARD 8: FORMAT (A4) 

If this card is blank and 

1 END was read at the end of the data deck a new 
data deck is read starting from Card 1 

2 STOP was read at the end of the data deck, the 
job is complete 

IE this card is KEPT the input will be repeated from Card 2 
and the same structure with different plotting parameters will 
be plotted. 
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SECTION 7 


EXAMPLE INPUT 



EXAMPLE INPUT 


The following pages contain the input data decks for a number 
of sample problems. The principal intent in providing sample in- 
put is to aid the user in better understanding the discussion in 
the main body .of the users' manual. Also provided with each prob- 
lem is a brief discussion and representative data so that bench- 
mark results are available for comparison. A detailed picture of 

i 

the sample problem idealizations showing node and element labeling 
has not been provided. It is suggested that the input plotting 
program SATELLITE be used to generate computer plots of the input 
data, if desired. 

Many sample problems that demonstrate the capability of the 
methods used in the PLANS program can be found in Refs. 1~3. 
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PPPPPPPPPPPPP 
ppppppppppppp p 
ppp pppp 

ppp ppp 

ppp ppp 

ppp ppp 

ppp pppp 

pppppppppppppp 
ppppppppppppp 
ppp 
ppp 
ppp 
ppp 
ppp 
.ppp 
ppp 


LLL 
LLL 
LU- 
LL. L 
LLL 
LLL 
LLL 
LLL 
LLL 
LLL 
LLL 
LLL 
LLL 
LLL 

llllllllllllll 

LLLLLLLLLLLLLL 


AA 

AAA A 
AA AA A A 
AAA AAA 
AAA AAA 

AAA AAA 

AAA AAA 

AAA AAA 


A AA AA AAAA AA AAAA A 
A AAAAAAAA AAAAAA A 
AAA AA A 
AAA AAA 
AAA AAA 
AAA AAA 
AAA AAA 
AAA AAA 


NNN 

NNNN 

NNNN 

NNNNN 

NNNNNN 

NNN NNN 

NNN NN N 


NNN 

NNN 

NNN 

NNN 

NNN 

NNN 

NNN 

NNN 

NNN 


NNN 

NNN 

NNN 


NNN 

NNN 

NNN 

NNN 

NNN 

NNN 

NNN 

NNN 

NNN 

NNN 


NNN NNN 
NNNNNN 
NNNNN 
NNNN 
NNNN 
NNN 


ssssssssss 
ssssssssssss 
ssss ssss 

ssss sss 

sss 
ssss 
ssssss 
ssssssssss 
sssssssssss 
ssssssss 
ssss 
sss 

sss ssss 

SSSS SSSS 

ssssss sssssss 
sssssssssss 


GGGGGG GGGGGGGGGGGGGGGG GGGGG GGG G GGGGGGGGGGGGGGGGGGGGGGGG GGGG GGGGGGGGGGGGGG GGGGGG GGGGGGGGGGGGGGGGGGGGG 

gggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggggg 

GG GGGGGGGGGGGGGG GGGGGG GGGGGG GGGGGGGGGGGGGGGGGGGGG 
GGGGGGGGGGGGGG GGGGG GGG GGGG GGGGG G GGG G GG GG G 
GGGGGGGGGGGGG GGGGGGGGGGGGGGGGGGGGG 
GG GGGGGG GGGGGGGGGGGGGGGG GGG 
GGGG GGGGGG GGGG GGGG GGGG 

, GGGGGGGGGGGGGGGGGGGGG 

V GGGGGGGGGGGGGGGGGGGGG 

nj GGGGGGGGGGGGGGGGGGGGG 

GGGGGGGGGGGGGGGG GGGG 

GGGGGGGGGGGGGGGG GG 

GGGGGGGGGGGGGGGG 

GG GGGGGG GGGGGG 

GG GGGGGG GGG6 

G G GGGGGGGG 

GG GGGGGG 

GGGGGG 

GGGG 

GG 


BEND MODULE 


ELASTIC-PLASTIC ANALYSIS 
PREPARED BY 

THE RESEARCH DEPARTMENT 
OF GRUMMAN AEROSPACE CORPORATION 
UNDER PARTIAL FUNDING FROM 
THE NASA LANGLEY RESEARCH CENTER 



Example Problem No. 1 
Program: BEND 

Problem Title: Classical Truss - Axial Force Element 

Comments: This problem uses the stringer element STGR in a 

simple 4 element truss. Symmetry conditions were used so that 
only half the truss was modeled. 

The accompanying figure shows the applied load versus the 
force in the vertical element (element 3). These result's were 
generated using a load step of AP = 848.5 lbs or 3774.3 N 
(which is 0.01 of the initial yield load). Agreement with 
results from an exact analysis was quite good. 
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Bend Sample Problem # 1 Classical Truss STRG Axial Force Elements 
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PROGRAM LISTING OF INKfT DATA CARDS 


e*b*oo*#*l»«t»«e*»ai^*D»*t»o*a«?oA(i«»»»a«4fiD»*»v0c«S«*«a0«*<to6ft«a . . • <n i>7«atk»aiaaR 

1234 567 890 1 234567890 12345678901 234 5 67890 1234 56789 01 2345678901 2345678901234667890 


I 

Ui 


BEND SAMPLE PROBLEM 41 CLASSICAL TR 

0 50 

1 ~4 




STKG 

1 

1 

2 



STRG 

2 

2 

4 



STRG 

3 

2 

3 



STRG 

SEND 

4 

1 

4 



0.0 


1 




.5 


2 

3 



1.0 


4 


BLANK 

CARD 

0.0 


1 

4 3 



.5 


2 


BLANK 

CARD 

0.0 


1 

-4 

BLANK 

CARD 

000000 


1 

3 



110000 


2 


BLANK 

CARD 

100000 


4 


BLANK 

CARD 

MSTG 1 .0 


EG 07 

1.0 E&00 

BLANK 

CARD 

1 -3 






MSTG 1.0 


£607 

0.5 FGOO 

BLANK 

CARD 

4 






SEND 

CONC 






4 0.5 


EG 00 

0.0 EGOO 

BLANK 

CARD 





BLANK 

CARD 

SEND 






1 ~4 




BLANK 

CARD 

-4 




BLANK 

CARD 

250000.0 


.Ol 



BLANK CARD 

STOP 


USS ST RG AXIAL FORCE ELEMENTS 


10.0 

EGOO 

1*0 

EG 05 

6.0 

EG 04 

10.0 

EGOO 

1*0 

EG 05 

t .0 

E&04 

o 

* 

o 

EGO 0 






t> * 1 s o . • o o 2 o o a *a**3*aaav b »4 m o a * » » 6 A • t> m * 6 * m 7 a o*« *8 

12345678901 234 5678 90 1234 667890 1234 5 67890 12345678901 2345 678901 2345678901234567890 



Example Problem No. 2 
Program: BEND 

Problem Title: NASA Frame with an End Load - BEAM ELEMENTS 

Comments:' This problem demonstrates the use of beam elements 

in an open space framework. Fifty seven beam elements of various 
cross section (L- section, T-section, rectangular section) con- 
necting 40 nodes are used in the model. Multipoint constraints 
were used to enforce symmetry conditions between the left and 
right hand side of the structure. Ideally plastic behavior was 
assumed with a yield stress of 50130 psi or 345.6 MPa. 

The figure shows displacement in the vertical and axial di- 
rections at the first bulkhead at the centerline (node 108) and 
at the midpoint between the first two bulkheads in the upper 
stringer (node 1001). The results indicate a plastic collapse 
at a load of P = 15000 lbs or 66.7 RN. 
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40 
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Bend Sample Problem #2 NASA Frame with End Load - 


KIPS 
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.P.RC6RAM .t_I STING OF 1 NPJT DATA CA 1 ? DS 


......... 1 2 J A 5 6 7 Q 

1234S67S9C123456 7B901234 5to78 9012345o78?:>1234 56 7S9D 1 2 34567890 1 2 3 45678 901 2 3456 7890 


BEND SAMPLE PROBLEM =2 h AS A FRAME WITH AN END LOAD - BEAM E_ £M ENTS 


-vj 

I 

00 


l C 1 

1CC1 102 

1 05 

2 C 1 2000 

108 3001 

208 

3 C 03 1 0C5 

2 06 

304 6002 

311 403 

4 I 1 

400 4CC4 

414 

80 02 

422 

8 OC 4 406 

BcAM 

1 

1 01 

1 02 

1 OB 


BEAM 

3 

1 06 

1 02 

1 08 


BEAM 

4 

108 

1 06 

1 02 


BEAM 

5 

201 

2 02 

208 


BEAM 

7 

206 

202 

20b 


BEAM 

8 

208 

2 06 

202 


BEAM 

1 0 

303 

3 04 

316 


BEAM 

14 

31 2 

3 04 

316 


BEAM 

16 

31 6 

312 

304 


BEAM 

18 

400 

4 08 

4 06 


BEAM 

22 

422 

4 08 

4 06 


BEAM 

24 

4 1 4 

422 

4 Ob 


BEAM 

25 

4 14 

4 06 

422 


bcA M 

27 

406 

4 00 

422 


b£AM 

.29 

1 02 

2 0 C 0 

1 Oo 


BEAM 

30 

2000 

2 02 

1 06 


BEAM 

31 

202 

20 02 

1 06 


BEAM 

32 

2002 

-3 04 

1 06 


BEAM 

33 

3 04 

2004 

106 


BEAM 

34 

2004 

4 08 

1 06 


BEAM 

35 

4000 

1 06 

1 02 


BEAM 

36 

2 06 

4000 

1 02 


BEAM 

37 

4002 

2 06 

102 


BEAM 

38 

312 

40 02 

1 02 


BEAM 

39 

4004 

3 12 

102 


BEAM 

40 

422 

4004 

1 02 


BEAM 

41 

6000 

1 08 

31 1 


BEAM 

42 

2 08 

60 00 

311 


BEAM 

43 

6002 

208 

31 1 


BEAM 

44 

316 

6002 

311 


BEAM 

45 

6004 

3 1 6 

311 


BEAM 

46 

4 14 

6004 

311 


BEAM 

47 

8000 

316 

311 


BEAM 

48 

80 02 

8000 

311 


BEAM 

49 

8004 

80 02 

31 1 


BEAM 

50 

4 06 

8004 

311 


BEAM 

53 

1 05 

1 01 

1 08 


BEAM 

54 

1 05 

1 08 

101 


BEAM 

57 

2 05 

201 

206 




»•••<£• 


• « • .3 • 


1234567d9012i4567«S/Ql 

2345o 76901 

2 345o7d90 


10 0 3 
ADO 2 


106 

2004 


202 

316 


6000 

3005 


205 

312 


303 

408 


4000 

6004 


2002 

8000 



PROGRAM LISTING -JT INPJT DATA CARDS 


1 3 . * . 4 . . . ......... 6 ... . .... . S 

12 34 36 7 3901 2 34567390 12 345678901 2345 07393 1234 567393 12 3456 78 901 23456789012 34567890 


•4 

i 

D 


BLANK CARD 


BcAM 


58 

2 05 

2 08 

BEAM 


64 

311 

3 03 

BEAM 


66 

31 1 

3 1 6 

8 LAM 


68 

403 

4 00 

BEAM 


72 

4 1 1 

4 03 

BEAM 


73 

41 1 

4 1 4 

BEAM 


75 

1 C01 

1 01 

BEAM 


76 

201 

1001 

BEAM 


77 

1 003 

2 01 

BCAM 


78 

3 03 

1 0 03 

BEAM 


79 

1005 

303 

BEAM 


80 

403 

1005 

BEAM 


81 

1 05 

3001 

tscAM 


82 

3001 

£ 05 

BEAM 


63 

2 05 

30 03 

BEAM 


84 

3003 

3 11 

BEAM 


85 

31 1 

3 0 05 

BEAM 


86 

3005 

4 1 1 

SEND 


56.875ECC 

101 

1 02 



39. S75E CC 

201 

202 



22.625E 00 

3 03 

3 04 



C.C E+OC 

403 

4 06 



4 6.3 75E CO 

2COO 

40 00 



31 . 25EC0 

2C 02 

40 02 



1 1 . 31 2 5EC0 

2 0 04 

40 04 


1 6 o 968 75E CC 

3000 




£ . 65625E CO 

80 04 




C.C EfOO 

1 08 

208 



C.C C+ CO 

6 0 04 

80 00 

1 

2 

• 3 2 3 1 2 5F. C 0 

1 02 

1 06 

-1 2 

. 3281 25EC0 

101 

1 OS 

1 

6 

. C31250ECC 

2 02 

206 

- 1 6 

. C31250E CO 

2 01 

205 

1 9 

• 

7890625EC0 

3 04 

3 12 

- 1 9 

• 

7690625600 

303 

3 1 1 

- 

24 . 71875E00 

403 

4 1 1 


24 . 71875E CO 

408 

422 

1 4 

• 

1 7963 75EC0 

2000 

40 00 

1 7. 

91 C 1 3625E CO 

2002 

40 02 

22. 

2 

£2 506256 CO 

2 C C4 

4 0 04 

- 1 4 

• 

1 796875E CO 

1C01 

3001 

- 1 7. 

91 Cl 5625E CO 

10 03 

30 03 


201 

315 

003 

4 06 
4 0b 
406 
105 
1 05 
1 05 
1 05 
1 05 
1 05 
101 
101 
101 
101 
1 01 
101 


1 05 

106 

105 


205 

2C6 

20 3 


31 1 

312 

315 


40S 

41 1 

414 

422 400 

1001 

30c : 

6000 


1003 

3003 

6002 


1005 

3005 

6004' 

800 2 


3 16 

400 

40 5 

414 6000 

80 02 

8004 

6002 



. , b. ...... . .7, ...... , . 8 

<i 34 567890 12 3456789012 34567890 


l 2 34 56 78901 2 345o7b 901 ; 
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PHCGGAM LISTING UF INPJT DATA CAkDS 


• • i ■ i >1 •• • • • g • • >2 t J • • • • • • • • >4 . • •• • • ■ . >3 .• . •••• • >6 i t o • •• • i •?•• •• •> • ■ ■ g 

1 2 34 5 6 7 6 901 2 3a5o 7890 1 2 3 4 56 7393 1 234567690 1234567390 12 34567390 1 2345678901 2 34567890 


BLANK CARD 


BLANK CARD 


BLANK CARD 


- 22. 25 J9C625E CC 

1 CCS 

30 Cb 


C.C E+CO 

406 



1 C . 1 Cl 5625E 00 

1 01 

1 02 


-10.1 Cl 5625E CO 

1 05 

1 06 

106 

1 3 . 1 55 31 25E 00 

201 

2 02 


- 13 . 1 9531 2 5S CO 

2 05 

2 06 

203 

16 . 32 C 31 2 5E CO 

3 03 

3 04 


- 16 .32031 25E 00 

311 

312 

316 

-3. 16C15625EC0 

8002 


■ 

2C.43 75E CC 

4 00 

403 

408 

2 0.43 75E 0 0 

411 

4 14 

422 

1 1 . 64843 75E CO 

2000 

1 001 


14 . 75781 25E CO 

20 02 

1 003 


1 3. 37890625E CO 

20 04 

1 005 


- 11 .64843751" CO 

4000 

3001 

6000 

- 1 4 . 75781 25E 00 

4C 02 

3003 

6002 

- 1 8. 3 7890625E 00 

4C 04 

30 0b 

60 04 


- 12. 24C23438E CC 8000 

- * . 0 8CC781 25E CO 8004 


1 01 0 1 C 


1 01 

2 01 

303 





1 01 0 1 C 


6000 

6002 

60 04 

8000 

8002 

800 4 

103 

1 Cl 0 1 C 


1 OS 

205 

311 





COOO CC 


400 

4 03 

408 

406 

411 

4 14 

422 

0000 CC 


1 02 

1 06 

202 

206 

304 

31 2 


CCOO CC 


1 001 

1003 

1005 

300 1 

30 0 3 

300 5 


10 2 

1 

1 01 

1 

1 .0 





1 06 

1 

1 05 

1 

1 .0 





102 

3 

101 

3 

1 .0 





i o e 

3 

1 05 

3 

1 

. . 0 




102 

5 

1 01 

5 

1 .0 





106 

5 

1 05 

5 

1 .0 





20 2 

1 

201 

1 

1 .0 





20 6 

1 

2 05 

1 

1 .0 





202 

3 

201 

3 

1 .0 





206 

3 

205 

3 

1 .0 





202 

5 

201 

5 

1 .0 





206 

5 

2 05 

5 

1 .0 





304 

1 

3 03 

1 

1 .0 





31 2 

1 

311 

1 

1 .0 





304 

3 

3 0 3 

3 

1 .0 





31 2 

3 

31 1 

3 

1 .0 






• •••••tlillll >lllbll«| l lll l Olll><>ll>7l> 6 

1 i34 56 7890 1 2 345676 90 1 2^>4 5673 90 12 34567390 1234 56 78 yO 12 34 067 84 Cl 2345678901 2 34 56 7 39 0 



PROGRAM LISTING U(- INPUT DATA CARDS 


...1 5 ........ .6 7 8 

1 234 567890123456789012345679901 2 34 5 673 93 1234 5673 90 12 34567 690 I 2 3 4567 S 9 01 2 34 56 789 C 




304 

5 

303 

5 

1 .0 








31 2 

5 

31 1 

5 

1 .0 








1 00 1 

l 

2000 

1 

1 . 0 








1003 

1 

2002 

1 

1 . 0 








100 5 

1 

20 04 

1 

1 . 0 








300 1 

1 

40 0 0 

1 

1 . 0 








300 3 

1 

4002 

l 

1 . 0 








3005 

1 

4004 

1 

1 . 0 








100 1 

2 

2000 

2 

, -1.0 








1003 

2 

2002 

2 

i -1,0 








100 5 

2 

2004 

2 

-1 . 0 








300 1 

2 

4000 

2 

-1 . c 








3003 

2 

4002 

2 

-1 . 0 








3005 

2 

4004 

2 

-1 . 0 








100 1 

3 

2000 

3 

1 . 0 








1003 

3 

2002 

3 

1 . 0 








100 5 

3 

2004 

3 

1 .0 








300 1 

3 

4 000 

3 

1 .0 








300 3 

3 

4002 

3 

1 . 0 








3005 

3 

4004 

3 

1 , 0 








100 1 

4 

2000 

• 4 

-1 .0 








1003 

4 

2002 

4 

-1.0 








1 00 5 

4 

2004 

4 

-1 .0 








300 1 

4 

40 00 

4 

-1 . 0 








300 3 

4 

4002 

4 

-1.0. 








3005 

4 

4004 

4 

-1 . 0 








100 1 

5 

2000 

5 

1 . 0 








1003 

5 

2002 

5 

1 .0 








1005 

5 

2004 

5 

1 .0 








300 1 

5 

4000 

5 

1 .0 








3003 

5 

4C02 

5 

1 .0 








, 300 5 

5 

4004 

5 

1 .0 








100 1 

6 

2000 

6 

-1 . 0 








1 003 

6 

2002 

6 

-1 .0 








1 00 5 

6 

2004 

6 

-1 . 0 








300 1 

6 

4000 

6 

-1 .0 








300 3 

6 

4002 

6 

-1 . 0 








3005 

6 

4004 

6 

-1 . 0 






BLANK 

CARD 

MdM 


1 1 . 2E0 6 


.333 JEOO 





5 Cl 30.E00 

BLANK 

CARD 

1 

3 

4 75 

7 

3 10 

14 

16 18 

22 

24 2 5 

27 29 -50 



5 3 

54 

57 56 

64 

66 63 

72 

7 3 75 

-86 

5 




SPEC 


2. OEOO 


. 1D6 7E00 


.5667EOO 


0 .OEOO 

.458 


......... 1 . . ...... .2 .3 . ....... .4 ... 5.... b . . ...... .7. ^ 8 

1 234 567 8901 2 3 4567890 1 2 3 4567893 12 34 5 67890 1 234 567890 1 2 34-567890 1 234 567 8 901 2 3456 789 0 



F&GGPAM LISTING OF 1 NP J T DATA CARDS 


.1 . .4 a 6 ......... 7, . . g 

1 2 34 56 78901 2 34567890 123.456 7890-1 234567890 1 234 567890 12345678901 234 5678901 2 34567890 


'•J 

I-* 

NJ 


BLANK CARD 
BLANK CARD 


0.0EC0 O.OEOO 

2.0EC0 l.OECO 

13 4 5 7 3 

54 57 58 64 66 

SPEC 6.0200 

0. OE 00 
6. OE CO 
25 , 27 

TSEC • . 4474E 0 0 


0 . 5 EOO 
O.OEOO 
1 . OEOO 


10 14 1 6 -IB 22 

■ 1 8 .OEOO 


tse: 


c. o 

i .-c 

4 1 V 46 


.90625 

2.09 


. 4474E00 


1 . 0 . 

4 7 -50 

L SEC 

. 93 75 
1 . C - 

29 -40-' 75 -86 

SEND 
CONC 


2.09 

.1211E00 .01099 
.9375 
1 ..0 


.0448331 .0729367 

.2775 

.04483 3 1 .0 7 29 36 7 

.2775 
.01099 
.0 625 




1 C 1 

1 3 

BLANK 

CARD 


t ‘ 



102 

t .. L 3 

BLANK 

CARD 

10 5 

- . i‘3 

BLANK 

CARD 

1 06 

-. 1-3 

BLANK 

CARD 

1 06 

- . 48 

BLANK 

CARD 



BLANK 

CARD 

SEND 




2 9 

.-50 

BLANK 

CARD 

-900 0 



5CC00. 0 


,01 


0 .0 

.0 9375 

0 .0 

,09375 
- .007037 
»0 62 5 


>4 6 8 72 73 S3 

O.OEOO 1 » 79,1 

7. 0293E-03 


7. 0293E-03 


1 .57675E- 04 


STOP 


..1 .......8 

12 3*+ 5678 90 123 45 o78 9 0123 *55 73 901234 5 67897123456789012 34 So 78 90 1234 56 7 89 012 34 66 7 89 0 



Example Problem No. 3 
Program: BEND 

Problem Title: 8 Element Square Plate 

Comments: This problem represents a simple statically deter- 

minate case that uses all combinations of the triangular membrane 
element, TRIM. These range from 3-node constant strain triangles, 
the 6-node linear strain triangle, and the two interface cases 
with 4 and 5 nodes. A distributed edge load (using the TRIA 
feature) is applied on one edge. 

The figure shows load versus longitudinal strain for a linear 
strain hardening material. 
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PRCGRAM LISTING O f= INPUT CATA CARDS 


„ . . 2 4 £ 6 ( • t t7i 

1 234567 B90 1234 56 7.8901 2 34 5678 90 1 234567890 1 234567890 l 2-345676901 2345678901 23456789 


'-J 

l 

(— 1 
Ln 


BLANK CARD 
BLANK CARD 
BLANK CARD 

BLANK CARD 
BLANK CARD 


BLANK CARD 


BEND SAMPLE PROBLEM =3 8 ELEMENT SQUARE PLATE TRIM ELEMENTS 
1 1 


lot -105 

- 202 

-205 

30 1 

303 

305 

404 

4C5 501 503 

£ 05 

TRIM 

30 1 

30 t 

0 

303 

0 

ec 2 

0 


TRI M 

40 1 

5 03 

0 

501 

0 

30 l 

0 


TRI M 

20 1 

303 

0 

30 l 

0 

10 1 

202 


TRI M 

402 

505 

0 

503 

0 

303 

404 


TRI M 

202 

3 05 

0 

303 

20 3 

103 

204 


TRI M 

302 

3 03 

0 

305 

405 

505 

404 


TRI M 

1 0 1 

101 

102 

103 

203 

303 

202 


TRIM 

SEND 

102 

1 03 

104 

l 05 

1 

205 

305 

204 


0 . 0 

E 00 

10 1 

30 i 

50 1 





5. 0 

E 00 

103 

303 

503 





1 0. 0 

E 0 0 

105 

305 

5 05 





o 

• 

o 

E OC 

101 

103 

105 





5.CCCCE 0C 

301 

30 3 

305 





1 C . 0 0 0E 0 C 

50 1 

503 

505 





0 . 0 

E 00 

10 1 

103 

1 05 





C.C0C0F 00 

301 

303 

305 





C.CCC0E OC 

501 

503 

505 





ccoooc 


101 







1 0000 c 


1 02 

-105 






1 1 coo c 


202 

-20 5 

303 

305 

404 

405 503 505 


CI000C 


301 

50 1 






vati i.o 


E 07 

1.0 


E 07 


0.0 2.0785 E 06 

1*0 

E-06 

l 

.0 

E-06 





3.6 

E 04 

3.6 


E 04 

3 .6 


E 04 2. 07 65 

E 04 

• 1 



0 .0 


1.0 


E 07 3.6 

E 04 3.6 

101 IC2 

20 1 

202 

30 1 

302 

401 

402 



thik i . 

0 








101 102 

20 I 

202 

30 1 

302 

40 1 

402 



SE ND 
TRI A 









501 

503 

3 01 

0 

.0 


_ 

1 .0 




0.0 



-1 

.0 



503 

505 

303 

c 

. 0 


— 

1 .0 




0.0 



-1 

.0 




. 3 


E 04 


• 0 « • « 4 • « tt 1 • • • « d • I • o2e D l • 4 0 4 » o3«44<l«l f • • 4 » O * • * 0 4 4 |5* « *• « • O 0 9 6« •4444t4«7o44Q4V44o8 

1 234-567 89 0 12 34 5678901 234667890 1234567890 1 23456 7890 t 234567 89 0 1 2345678901 234567690 


ac o 



PROGRAM LISTING OF INPUT DATA CARDS 


4 ..5.... 6o... 7... 8 

1-2 34567 890 1234 56 789012 34567890 1 234567890 1 2-34567890 1 23456789012 3456 78901 234567890 


BLANK CARD 

BLANK CARD 
BLANK CARD 


SET NO 
-1 000 

-1000 

5CCC0.0 


STOP 


.0 l 


l 


On 


. . , « « a . . . 1 . . . 2. ........ ....... ..5. .6. .8 

1234567890 1234 56 78901234 567890 1 234567890 12345678901 234567 89012 345678901 234567890 



Example Problem No. 4 
Program : BEND 

Problem Title: Collapse of a Uniformly Loaded Circular Plate 

Comments: This problem demonstrates the use of the plate ele- 

ment TRIP. One quadrant of a uniformly loaded, simply supported 
circular plate was modeled consisting of 50 elements, 36 nodes, 
and 151 degrees of freedom. The multipoint constraint feature 
was used to transform the edge rotations to their normal’ and tan- 
gential components in order to apply the simple support boundary 
conditions (0^ =0). To be noted also is that the order of the 
element input was set so that the local directions of the curva- 
tures at each node (defined by the first element containing that 
node) line up with the plate boundaries and lines of symmetry. 

This was done to enforce boundary/ symmetry conditions for the 
higher order degrees of freedom (see Appendix IV, Section 2). 

A complete discussion of the plastic analysis of this plate 
as well as a number of other plate configurations is given in 
Ref. 2. The figure shown here is for the load versus central de- 
flection up to the collapse load. 
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LOAD VERSUS CENTRAL DEFLECTION 


Bend Sample Problem #4 Uniformly Loaded Simply Supported Circular Plate 
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program listing of input data cards 


. 1 ........ .2 A 5 .. 6 8 

1 2345678901 2 34 567 8901 2 34 56 7890 1 23 4 56 789 0 1 2345678901 2345678901 2345678901 234567890 


PEND £A^LE PRCPLPN' =4 UNIFORMLY L.QAC5C CIRCULAR PLATE 


1 

1 0 

1 

1 

50 


1 

TRI <= 

-36 

1 

2 

1 

3 

TP I a 


2 

1 

4 

3 

TRI P 


3 

5 

2 

6 

TRI P 


4 

2 

3 

6 

TP I P 


5 

7 

6 

3 

TRI P 


6 

8 

7 

3 

TP I P 


7 

3 

4 

8 

TRI P 


8 

4 

9 

8 

TRI P 


9 

10 

5 

1 1 

TP I P 


1 0 

5 

6 

1 1 

TRI P 


1 1 

1 2 

1 l 

6 

TRI P 


l 2 

6 

7 

12 

TRIP 


1 3 

13 

1 2 

7 

TRI P 


l 4 

14 

1 3 

7 

TRI P 


15 

7 

8 

14 

TRI P 


i e 

15 

1 4 

8 

TRI P 


1 7 

3 

9 

15 

TPI P 


i e 

9 

1 6 

1 5 

TRI P 


t 9 

17 

10 

18 

TRI P 


2 C 

1 0 

1 1 

18 

TRI P 


2 1 

19 

1 8 

1 1 

TF I P 


22 

1 1 

1 2 

19 

TRI P 


23 

20 

19 

1 2 

TtRl P 


24 

1 2 

1 3 

20 

TRI P 


25 

21 

20 

■ 13 

TPI P 


26 

22 

2 l 

13 

TRI P 


27 

13 

1 4 

22 

TRI P 


28 

23 

22 

14 

TRI o 


29 

14 

15 

23 

TRI P 


30 

24 

23 

1 5 

TPI P 


3 1 

15 

16 

24 

TR I n 


32 

16 

25 

24 

TRI P 


33 

26 

1 7 

27 

TPI P 


34 

1 7 

18 

27 

TPI P 


35 

28 

27 

18 

TPI P 


3 6 

18 

19 

28 

TRI P 


37 

29 

28 

19 

TRIP 


3 8 

19 

20 

29 

TRI P 


39 

30 

29 

20 

TRI P 


4 C 

?0 

2 1 

30 

TRI P 


4 1 

31 

30 

21 


...» 1 ...2. ...... ..3... 4...o 5 6. ...... ..7 

1 234 5 67 89 0 12 34 56 78901? 34 567890 12 34 56 7890123456 78 901 2 34567 8 9 0 1 2 3456 789 01 23456 789 0 
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PROGRAM LISTING OF INPUT DATA CAROS 

1... 2 3. 4 8 

123456789012345678901 2345678901234567890 1 22456 7890 1 23456 7 890 1 2 345678901 234567890 


BLANK CARD 


TRI P 

4 2 

32 

3 1 

21 


TPI P 

43 

21 

2 2 

32 


TRI P 

44 

33 

32 

22 


TRI P 

45 

2? 

23 

33 


TRI P 

46 

34 

33 

23 


TRI P 

47 

23 

24 

34 


TRI P 

50 

25 

36 

35 


TRI P 

48 

35 

34 

24 


TRI P 

49 

24 

25 

35 


SEND 






0. 0 


1 

2 

5 

1 0 

1.41422 


3 

6 

1 1 

18 

2.0 


4 




2.82844 


7 

1 2 

19 

28 

3. 74162 


a 




4.0 


9 




4.24266 


13 

20 

29 


5.29146 


14 




5. P3C92 


15 




6. 0 


16 




5.65668 


21 

30 



6.78232 


22 




7.46326 


23 




7. 874 


24 


. 


8. 0 


25 




7. C71 1 


31 




8.24621 


32 




9.05538 


33 




9. 59166 


34 




9. 89949 


35 




10.0 


36 




0.0 


1 

4 

9 

16 

2 . 0 


2 




1 . A 1 422 


3 

8 

15 

24 

4 . 0 


5 




3.74162 


6 




2. 82844 


7 

1 4 

23 

34 

6.0 


10 




5. 83 C92 


1 1 




5.29146 


1 2 




4.24266 


13 

22 

23 


8. 0 


17 




7.8 74 


18 





17 

27 


25 

35 


26 


36 


. . . . 1.. 2 .3... 4 ....5. ...... ..6... ....... .8 

123 45 6 789012 34 56 78901234 56789 0 I 234 56 789 0 1 224 56 78 90 1 234567 690 1 2 345678901 23456 7890 
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program listing OF input DATA CAROS 


BLANK CARD 
BLANK CARD 


BLANK CARD 


BLANK CARD 


1 . . 2. ...... . .3 4 5. ...... ..6.... 7 3 

12 345 67 89012 34 56 78901 2.34567390 123456 78 90 1 224567890 1 234567890 12 3456 78 901 23456 7890 


7. 43328 


19 



6. 78232 


20 



5.65688 


21 

3 2 


10.0 


26 



9 . 8994 9 


27 



9. 59166 


23 



9. 0553 8 


29 



8,2462 1 


3C 



7. C7 1 1 


31 



0. 0 


1 

-36 


C01 1 1 <-0001 1 1 

1 

-36 


COl C0CCCC1 1C 

1 



CC11 CCCCC1 10 

2 

5 

10 17 

CC1 Cl CC CC 1 1 c 

4 

9 

16 25 

C OC 01 COCOOOl 

27 

-35 


COCCI ccoocoo 

26 

36 


26 

4 

26 

5 

-1 .0 

26 

« 

2 6 

5 

0.0 

2. 7 

4 

27 

5 

-.98994904 

2 7 

5 

27 

5 

.14 142 199 

28 

4 

28 

5 

- .9 59 16595 

2 8 

5 

28 

5 

. 2828439 7 

29 

4 

29 

5 

-.90553799 

29 

5 

?9 

5 

.42426596 

3 0 

4 

30 

5 

-.82462 10 1 

30 

5 

30 

5 

.565638 C4 

31 

4 

71 

5 

-.7C711C C2 

31 

5 

31 

5 

.707 1 1 CO 2 

32 

4 

32 

5 

-.56568804 

32 

5 

32 

5 

.82462 10 1 

3 3 

4 

37 

5 

- .42426556 

33 

5 

33 

5 

.90553799 

3 4 

4 

34 

5 

-.28284397 

34 

t) 

14 

5 

.95916555 

36 

4 

35 

5 

- . 1 41 42 199 

3 5 

5 

35 

5 

.98994504 

36 

4 

3 6 

5 

0 .0 

36 

5 

36 

5 

1 .0 

VAT1 

1 . 0 

K + 07 

1 . 

C F+0 7 O.C 4 • 03 22 5 3E + 06 .24 

C. 0 


0 

.0 



1 • * O • 4 

o 

» • • s • *1 • 



121456789 01234 56 78901 2 34 5 6 78<>0 1 23 4 5678 9 0 1 2 34 5 6 789 C 1 234 56 7 £9 3 1 2 34567 89 0 1 23456 789 0 



7-22 


BLANK CARD 
BLANK CARD 

BLANK CARD 


PFGGPAH L'IST ING CF INPUT DATA CAPOS 

1 • > i • . • • 3 . > « < i > . . • A < 5 1 •••••••• 6 ••<><••(. 7 • • ( 6 

1 234567890123 4"5f> 7 8 901 2 34 56789 0 123456 7890 1234567890 1 2 345 57 P90 12 3456789 01 2 3456789 0 


36CCC .0 
C. 6 

1 -50 

THI K 1.0 

! 1 -5*0 

NLRS 

1 0 1 

SF ND 

SURF 0.0 

C. C 
C. 0 

1 -5 0 

SEND 

1 2 . 


36000 .0 ‘ 

0.0 


0 .0 


0 .0 

0 .0 - 


3600C. 0 
1.0 E+C7 


- 1.0 

- 1.0 


- 1 .0 


2C764.6 

36000.0 


36000.0 


1 -36 

630.0 


.0 1 


STOP 


1 2. 3. ....... 6... 7 8 

12 34 5 67 89 0 12 54 6678 901 ? 34 567890 1 2 3 4 56 7 89C 1 234St:7 8 9C l 234 567 6 9 0 1 2 345678 90 1 23456 7890 



Example Problem No. 5 
Program: BEND 

Problem Title: Thermoplastic Analysis of a Circular Disk 

Comments: This problem demonstrates the thermal stress capa- 

bility in conjunction with the triangular membrane element. To do 
this, an elastic plastic analysis of an annular disk subjected to 
a steady state axisymmetric radial temperature distribution was 
performed. Since the response is axisymmetric, only a slice of 

i 

the disk was modeled with the multipoint constraint feature being 
used to enforce tangential and normal boundary conditions on one 
edge. Forty constant strain triangular elements connecting 32 
nodes were used in the idealization. The steady state temperature 
distribution is applied at .nodes according to the relation 

T =_T a _- (T b - T a ) £n(r/a)/£n(b/a) 

where T^ and are the inner and outer radius temperature, 

respectively. 

Results shown in the figure are for dimensionless circumfer- 
ential stress versus the radial coordinate. Comparison between 
these results and Ref. 4 are good. 
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PROGRAM LISTING OF INPUT DATA CARDS 


• S«c0t>o29eo«»*»»e3«a»»**»«*4*«t«aii*a*5oft0d«o»o*6**«f>o**«e7ft«»*oa«ao8 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 


1 

-3 

7 

-9 

13 

-15 

1 9 

-21 

49 

-51 

55 

-57 

61 

63 



TRIM 


0 1 

01 

0 

03 

0 

02 

TRIM 


02 

03 

0 

09 

0 

02 

TRIM 


03 

09 

0 

07 

0 

02 

TRIM 


04 

07 

o 

01 

0 

02 

TRIM 


05 

07 

0 

09 

0 

08 

TRIM 


06 

09 

0 

15 

0 

08 

trim 


07 

15 

0 

13 

0 

08 

TRIM 


08 

13 

0 

07 

0 

08 

TRIM 


09 

1 3 

0 

15 

0 

14 

TRIM 


10 

15 

0 

21 

0 

14 

TRIM ' 


1 I 

21 

0 

19 

0 

14 

TRIM 


12 

19 

0 

13 

0 

14 

TRIM 


13 

19 

o' 

21 

0 

20 

TRIM 


14 

21 

0 

27 

0 

20 

TRIM 


1 5 

27 

0 

25 

0 

20 

TRIM 


16 

25 

0 

19 

0 

20 

TRIM 


1 7 

25 

0 

27 

0 

26 

TRIM 


18 

27 

0 

33 

0 

26 

TRIM 


19 

33 

0 

31 

0 

26 

TRIM 


20 

3 1 

0 

25 

0 

26 

TRIM 


23 

31 

0 

33 

0 

32 

TRIM 


22 

33 

0 

39 

0 

32 

TRIM 


23 

39 

0 

37 

0 

32 

TRIM 


24 

37 

0 

31 

0 

32 

TRIM 


25 

37 

0 

39 

0 

38 

TRIM 


26 

39 

0 

45 

0 

38 

TRIM 


2? 

45 

0 

43 

0 

38 

TRIM 


28 

43 

0 

37 

0 

38 

TRIM 


29 

43 

0 

45 

0 

44 

TRIM 


30 

45 

0 

51 

0 

44 

TRIM 


31 

51 

0 

49 

0 

44 

TRIM 


32 

49 

0 

43 

0 

44 

TRIM 


33 

49 

0 

51 

0 

50 

TRIM 


34 

51 

0 

67 

0 

50 

TRIM 


35 

57 

0 

55 

0 

50 

TRIM 


36 

55 

0 

49 

0 

50 

TRIM 


37 

55 

0 

57 

0 

56 

TRIM 


38 

57 

0 

63 

0 

56 

TRIM 


39 

63 

0 

61 

0 

56 

TRIM 


4 0 

61 

0 

55 

0 

56 

* «* o O « 0 • 

v a 1 © 

O t> <* * «> O 

» • 2» 


• »3o 


..4 


25 -27 


31 


-33 


37 -39 


1 234 o67890 1 234-55789 01 234 567890 1 234 567890 1 234-557 89 012345 678901 234567 890 


4 3 -4 5 


X* © 
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PROGRAM LISTING OF INPUT DATA CARDS 

1 2. . 7.«. 

12345678901 2395678901234 5678 901234 5 678901 23456789 01234 5 678901 2345678901 234567890 


BLANK CARD 


SEND 


1.0 

E 

00 

3 

6 

1.1 

E 

00 

9 

12 

1.2 

E 

00 

15 

18 

1.3 

E 

00 

21 

24 

1 .4 

E 

00 

27 

30 

1.5 

E 

00 

33 

36 

1.6 

E 

00 

39 

42 

1.7 

E 

00 

45 

48 

1.8 

E 

oc 

51 

54 

1.9 

E 

00 

57 

60 

2.0 

E 

00 

63 

66 

•99452189E 

00 

1 

4 

1.0939741 

E 

oo 

7 

10 

1. 1934 262 

E 

00 

S3 

16 

1.2928784 

E 

00 

19 

22. 

1.3923306 

E 

00 

25 

28 

1.491 7828 

E 

00 

31 

34 

1.5912350 

E 

00 

37 

40 

1.6814254 

E 

00 

43 

46 

1 . 790 1 394 

E 

00 

49 

52 

‘ 1.8895916 

E 

00 

55 

58 

1.9890437 

E 

00 

61 

64 

1.05 

E 

00 

2 

5 

1 . 15 

E 

00 

' 8 

1 1 

1 .25 

E 

00 

14 

17 

1 .35 

t 

00 

20 

23 

1.45 

E 

00 

26 

29 

1.55 

E 

00 

32 

35 

1 .65 

E 

00 

38 

4 1 

1 .75 

E 

00 

44 

47 

1.85 

E 

00 

50 

53 

1 .95 

E 

00 

56 

59 

0.0 

E 

00 

3 

-3 

.0523359 

E 

oo 

2 

5 

.0575695 

E 

oc 

8 

1 1 

.0628031 

E 

00 

14 

17 

.0680367 

E 

00 

20 

23 

.0732703 

E 

00 

26 

29 

. 0785039 

E 

00 

32 

35 

.0837375 

E 

00 

38 

4 1 

-088971 1 

E 

00 

44 

47 

.0942047 

E 

00 

50 

S3 


1 ... .3. . 6. ....••..7. ........8 

1 23 4 567890 1 2345t>789 01 234567890 1234 567890 1 23 4567 89 012345 678901 234567 890 1 234 567890 
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PROGRAM LISTING OP INPUT DATA CARDS 


• #9 » • tt B • • 1 Q • . 9CBBBb2b B B « B B6B»3bBBBB B B t B A« B B B B BIBb5bBBB« BB B b6b B B 9 9 9 9 Bb7bBBB4BB 0*3 

12345678901 234567890123456789012345676901 234567890 1 23456789012345678901234567890 



.0994383 E 

00 

56 

59 







.1046712 E 

00 

62 

65 







.1045284 E 

00 

1 

4 







.1149813 E 

00 

7 

SO 







. 1 25434 1 E 

00 

13 

16 







.1358070 E 

00 

19 

22 







.1463398 E 

00 

25 

28 







.1567926 E 

00 

31 

34 







.1672455 E 

00 

37 

40 







.1776983 E 

00 

43 

46 







.1881512 E 

00 

49 

52 







.1986040 E 

00 

55 

58 






BLANK CARD 

.2090569 E 

00 

6 1 

64 







0.0 


1 

r-63 






BLANK CARD 

1 10000000000 


1 

-3 

7 -9 

13 -15 19 

-21 

25 -27 31 -33 

37 




-39 

43 

-45 49 

-51 55 -5 7 

61 

63 


BLANK CARD 

100000000000 


1 

-6 

-61 3 

-6 -63 




1 1 


l 

1 

.9945218 

1 

2 

-.1045284 



I ' 2 


1 

1 

• 1045284 

1 

2 

.9945218 



7 1 


7 

1 

.9945218 

7 

2 

-.1045284 



7 2 


7 

1 

• 1045284 

7 

2 

.9945218 



S3 1 


13 

1 

.9945218 

13 

2 

-.1045284 



13 2 


1 3 

2 

. 1045284 

13 

2 

.994521 8 



19 1 


19 

1 

.9945218 

19 

2 

1045284 



19 2 


19 

I 

.1045284 

1 9 

2 

.9945218 



25 1 


25 

1 

.9945258 

25 

2 

-.1045284 



25 2 


25 

1 

.1045284 

25 

2 

.9945218 



31 1 


31 

S 

.9945218 

31 

2 

1045284 



31 2 


31 

1 

. 1045284 

3 1 

2 

.9945218 



37 1 


37 

1 

.9945218 

37 

2 

-.1045284 



37 2 


37 

1 

.3045284 

37 

2 

.9945218 



43 1 


43 

1 

.9945218 

43 

2 

-.1045264 



43 2 


43 

1 

. 1045284 

43 

2 

.9945218 



49 1 


49 

1 

.9945218 

49 

2 

— » 1045284 



49 2 


49 

l 

. 1045284 

49 

2 

.9945218 



55 1 


55 

5 

.9945218 

55 

2 

-.1045284 



55 2 


55 

3 

.1045284 

- 55 

2 

.9945218 



61 2 


6 1 

1 

.9945218 

61 

2 

-.1045284 



61 2 


6 1 

1 

. 1045284 

6 J 

2 

.9945218 


BLANK CARD 











MAT! 1.0 

E 

07 

1 

.0 E 07 

0.0 F 00 

3.846153 E 06 .30000E 

00 



» o » # 

» « 2* 

• « O O « 


• «o4« ® • m o o o • » 5 c 

0*000 


» 08 


1 234-667890 1 234 5670901 234 567890 1 234567890 1 234567890 1 2345 67890 1 234567 6901 234567890 



8 Z-L 


PROOF A V LISTING OF INPUT PAT A CARDS 


2 .*.. 5. ...... ..6 8 
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Example Problem No. 6 
Program : BEND 

Problem Title: Elastic Plastic Analysis of an SAE Keyhole Specimen 

Comments: The accompanying figures show an SAE keyhole speci- 

men, the finite element idealization of half of the plate, and 
results for the strain at the notch versus load. The idealization 
uses 92 constant strain elements connecting 64 nodes that leads 
to 120 degrees of freedom. 
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Bend Sample Problem #6 Load versus Notch Strain £>r SAE Keyhole Specimen 
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LOAD AMPLITUDE 
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16 r 



e a NOTCH STRAIN AMPLITUDE 

Bend Sample Problem #6 (Continued) 
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ELASTIC-PLASTIC ANALYSIS 


PREPARED BY 
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THE NASA LANGLEY RESEARCH CENTER 



Example Problem No. 7 
Program: REVBY 

Problem Title: Plastic Analysis of a Ring Stiffened Spherical Shell 

Comments: This problem demonstrates the use of the shell and 

ring element. The structure consists of a uniformly loaded spheri- 
cal shell with a stiffened circular hole at the shell apex. The 
problem is idealized using 26 axisymmetric shell elements and one 
ring element connecting 27 nodes. 

The figure shows the normal displacement versus applied 
pressure at the ring shell interface and also results for an un- 
stiffened shell. Sudden collapse occurs for the stiffened shell 
at qa^/Et^ — 15000. 

Additional results for this problem are in Ref. 3. 
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LOAD-DEFLECTION CURVES FOR RING-STIFFENED 
SPHERICAL SHELL UNDER EXTERNAL PRESSURE 



(a) NORMAL DISPLACEMENT AT THE HOLE BOUNDARY 


REVBY Sample Problem #1 Ring Stiffened Spherical Shell 
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PROGRAM LISTING OF INPJt DATA CARDS 


■-J 

1 

■P- 

L5 


BLANK 

BLANK 

BLANK 


BLANK 


'BLANK 


1234 367 89 01 2 3 A 56 78901 2 34 567 3 9 0 1 2 3 4 55 73 90 l 2345678 90 1 234567890 12 345678901234567890 



2. 4364647E-02 
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27 
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Ot CCCG0 0 05 0 0 

0.30 000 005 05 


1 . CC0C0C05-05 

1 . 00 00 0 0 05 — 0 5 
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*o.....a«lfl#***»».e2oo*.o©..*3*...«a.*.4©..o..o*a5..o.o..«.6«*.<>«..o*7..ooo..Q.8 

1234 5 67 8901234 56 78901 23456739012345573901 234567890 1 234567890 1 2 3456 7890 1 23456 7890 



PROG PA M LISTING OF INPJT OATA CAP CS 
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SEND 
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STOP 


• o o © • oMfi 1 • o • * 69i«i9 oe6*»*«*»ooQ7(iooe«09eo8 

12 34 567 89012 34 56 78901 2 34 567 3 901 2 3 455 739 0 1234567 390 1 234 567 890 12 345678901234567890 



Example Problem No. 8 
Program : KEVBY 

Problem Title: Thick Tube Under Internal Pressure 

Comments: This problem demonstrates the use of the axisym- 

metric revolved triangular element, TRIR. Five bays of four 
triangles each (20) are used to idealize thick tube under in- 
ternal pressure. 

A series of results are shown in the figure and compared with 
a solution in Ref. 5. Additional discussion of this problem is in 
Ref. 3. 


7-45 




PRESSURE p VERSUS RADIAL DISTRIBUTION OF RADIAL 

DISPLACEMENT u (b) AND RADIUS STRESS 

OF ELASTIC-PLASTIC BOUNDARY, p 



DISTRIBUTION OF CIRCUMFERENTIAL STRESS 

REVBY Sample Problem #2 Thick Cylinder Under Internal Pressure 
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Example Problem No. 9 
Program : REVBY 

Problem Title; Thermoplastic Analysis of a Circular Disk 

Comments: This problem is the same as sample problem 5 using 

the BEND program. Here 10 shell elements are used connecting 
11 nodes leading to 22 degrees of freedom. Comparison made with 
the results of Ref. 4 indicate excellent agreement. 
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PkOGRAM LISTING OF INPUT DATA CAROS 


A»««*ao*ol*oo*«ot)*o^a«o«»o««#3a4*»*o oft*4*o*a*o«a»5®a*#a*ooa6**»»B»«»6/oo*6**t**b 
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BLANK CARD 


BLANK CARD 
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REVbY SAMPLE PRU&lEM ft 3 THERMOPLASTIC ANALYSIS OF CIRCULAR DISK 

0 10 

1 -1 1 


SHEL 

1 

1 

2 

SHEL 

2 

2 

3 

SHEL 

3 

3 

4 

SHEL 

4 

4 

5 

SHEL 

5 

5 

6 

SHEL 

6 

6 

7 

SHEL 

7 

7 

8 

SHEL 

b 

8 

9 

SHEL 

9 

9 

10 

SHEL 

SEND 

10 

10 

11 

1 . OOOOOOOE 

00 

1 


1 olOOOOOOt 

00 

2 


1 .2000000E 

oo 

3 


1 ® 3000C0CE 

GO- 

4 


1 .4000000E 

GO 

5 


i .sooooooe 

oo 

6 


1 obOOOOOOt 

00 

7 


1 .7000000E 

00 

b 


1 .6000000E 

00 

y 


1 .9GCCOOOE 

00 

10 


2 ® OOOOOOOE 

00 

1 1 


O.QOOOOOOfc 

00 

1 

-11 

0 ® OOOOOOOE 

00 

1 

-11 


1 1 


lOOlii* 


1 -1 1 


10 .OOOOOOOE 

06 

lO. OOOOOOOE 
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PROGRAM LISTING Of= INPUT DATA CARDS 


• SoilIllMBsXSXIIlAottlKMSlIKMtMblaMdt.tTltdttOtb 
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12o 4 56 78961 2 34 56789 Cl 234 567890 1234 567690 1 2 34 fat. 78VG 1 2 34567 890 1 234 567890 1 234 56 7890 



program Lifting up inpjt data ca^dj 
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Example Problem No. 10 
Program : REVBY 

Problem Title: Thermoplastic Analysis of a Circular Disk 

Comments: Same as sample problem 5 (BEND) and 8 (REVBY). Use 

here is made of 40 revolved triangular element connceting 32 nodes 
leading to 53 degrees of freedom. Results shown in the figure are 
for average stresses at the intersection of the diagonals of each 
bay. Again comparison with the results from Ref. 4 are good. 
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FINITE ELEMENT IDEALIZATION USING REVOLVE TRIANGULAR ELEMENTS 



REVBY Sample Problem #4 Thermoplastic Analysis of a Circular Disk 
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Example Problem No. 11 

Program: HEX 

Problem Title: Bar Subjected to an End Displacement 

Comments: This is a simple statically determinate problem to 

demonstrate the use of the HEX element. The automatic mesh genera 

tion feature is used (MSGN) to specify a 1x6x4 mesh of 
24 elements. An applied edge displacement is imposed on one end. 

The figure shows the axial stress versus applied end displace 

ment. 
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PROGRAM LISTING OF INPUT DATA CARDS 


7 . ....... .8 

1 2345O7890 1 23456789 0123* 567890 123456 7890 123456 7890 12 345678901 234S678901 2 34 56 7890 


HEX SAMPLE PROBLEM ft 1 LONG BAR SUBJECTED TO AN APPLIED END DISPLACEMENT 

0 5 115 

1 -70 

MSGN 

16 4 11 

SEND 


0. OOOOOOOE 00 

1 

3 

5 

7 

9 

9 

11 

13 

15 

17 

19 

21 

23 

0. OOOOOOOE 00 

25 

27 

29 

31 

33 

35 

37 

39 

41 

43 

45 

47 

49 

0. OOOOOOOE 00 

51 

53 

55 

57 

59 

61 

63 

65 

67 

69 




2.5000000E-01 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

22 

24 

26 

2 . 50 00000 E—0 1 

28 

30 

32 

34 

36 

38 

40 

42 

44 

46 

48 

50 

52 

2 . 50 OOOOOE— 0 1 

54 

56 

58 

60 

62 

64 

66 

68 

70 






BLANK CARD 

0. OOOOOOOE 00 1 -10 

5.00 OOQOOE— 0 1 11 -20 

8. 7500000E— 01 21 -30 

1 . 25 OOOOOE 00 31 -40 

1.5000000E 00 41 -50 

1 a 75 OOOOOE 00 51 -60 

2.00 OOOOOE 00 61 -70 

BLANK CARD 


-5 . OO'OODOOE— 0 1 
—5. OO OOOOOE— 0 1 

1 

62 

2 

11 

12 

21 

22 

31 

32 

41 

42 

51 

52 

61 

—2 .50 OOOOOE— 01 
—2 • 50 OOOOOE —0 1 

3 

64 

4 

13 

14 

23 

24 

33 

34 

43 

44 

53 

54 

63 

0. OOOOOOOE 00 
0. OOOOOOOE 00 

5 

66 

6 

15 

16 

25 

26 

35 

36 

45 

46 

55 

56 

65 

2.50 OOOOOE— 01 
2 . 50 OOOOOE— 0 1 

7 

68 

8 

17 

18 

27 

28 

37 

38 

47 

48 

57 

58 

67 

5. 00 OOOOOE— 0 1 

9 

10 

19 

20 

29 

30 

39 

40 

49 

50 

59 

60 

69 


5.00 OOOOOE— 0 1 70 

BLANK CARD 

001 1379 

101 2 4 8 10 

000 5 

100 6 

121 61 -70 

61 2 70 2 1.0 

62 2 70 2 1.0 

63 2 70 2 1.0 

64 2 70 2 1.0 

65 2 70 2 1.0 

66 2 70 2 1.0 


....... a.l*. ...... .2. e. .o. *».3»o. . sc. ..4.s.ao. .o.5.»«. ..... 6.. a. .e..o7. ...... ..8 

1 2345678901 2345678901234567890 1234567890 1 234567890 1 234567890 12345678901 234567890 


BLANK CARD 
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PROGRAM LISTING OF INPUT DATA CARDS 


I 2 34 56 7690 12345 6789 01 234 567890 I 23456 7890 1 234567890 1 234 567890 1 234567890 1 234 567890 


67 2 70 2 1.0 

66 2 70 2 1.0 

69 2 70 2 1.0 

70 2 70 2 .003114 

BLANK CARD 

MAT I 3.0 E 07 3.0 E 07 3.0 E 07 .33333E 00 .33333E 00 


» 333333E 

00 

1.125 E 

07 

1 « 125 

E 

07 

1.125 e 

07 



0.0 E 

00 

0.0 E 

00 

0.0 

E 

00 





12. E- 

-06 

12. E- 

06 

12. 

E 

-06 





3.6 E 

04 

3.6 F 

04 

3.6 

E 

04 

2.0784E 

04 

2.0784E 

04 

2.0784 E 

04 

10. t 

00 

4.0 

E 

04 

3.6 E 

04 

3.0 E 

07 


1 -24 

SEND 

SEND 

-100 

BLANK CARD 

-J 00 

BLANK CARD 

1.5 .01 


BLANK CARD 


STOP 


.0.o.....l.o©o.«...2.o.....*.3......o.»4.*«.ooo.o5...i*....6...oo®...7».<>.ft.©.o6 

1 2 3 **56 7 6 90 1 23456 789 0 I 234 567890 1 284 567 690 1 234 56 78901 2 34 567890 ] 234 567S90 1 234 56 78 90 



Example Problem No. 12 
Program: HEX 

Problem Title: Thermoplastic Analysis of a Circular Disk 

Comments: This problem is the same as problems 5 (BEND), 

8 , and 9 (REVBY ) . Here the idealization of a slice of the disk 
is used with 10 hexahedra elements with midside nodes in the 
radial direction. Thirty Lobatto stress points are used in each 
direction, 10 in the radial, 3 in the circumference, and 1 in 
the thickness direction. 
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FINITE ELEMENT MODEL OF A SECTOR OF THE DISK 
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Example Problem No. 13 
Program: HEX 

f 

Problem Title: Collapse of a Simply Supported Uniformly Loaded Beam 

Comments: Ten 20 node hexahedra elements are used to model 

half of a simply supported beam. A 1x1x8' array of Lobatto 
points are used (8 through the thickness) to determine stresses 
within each element: Lobatto points were chosen in order to have 
stress points at the surface. Eight points were taken through the 
thickness in order to accurately define an elastic plastic boundary 
through the thickness. 

The figure shows the central deflection versus load for an 

ideally plastic material. Results are in good agreement with 

/ 

Ref. 5, 
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FINITE ELEMENT MODEL OF HALF OF THE BEAM 



.6 .7 .8 .9 1.0 P ~~2o '■? 

o 

LOAD VS. NON-DIMENSIONAL CENTER DEFLECTION (W Q CORRESPONDS TO 

DISPLACEMENT AT INITIAL YIELD) OF A SIMPLY-SUPPORTED BEAM. HEX 
RESULTS CORRESPOND TO AN1X1X8 ARRAY OF LOBATTO POINTS. 

HEX Sample Problem #3 Collapse of a Uniformly Loaded Beam 
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Example Problem No. 14 
Program: OPLANE~GM 

Problem Title: Uniformly Loaded Restrained Beam 

Comments: Ten beam elements are used to model half a simply 

supported restrained beam. A full tangent modulus method is used 
(PTAN = 0.0). 

The figure shows the central deflection and internal axial 
force versus the total load. 
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